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Abstract. We develop the kinetic theory of point vortices in two-dimensional hydrodynamics and illustrate 
the main results of the theory with numerical simulations. We first consider the evolution of the system 
"as a whole" and show that the evolution of the vorticity profile is due to resonances between different 
orbits of the point vortices. The evolution stops when the profile of angular velocity becomes monotonic 
even if the system has not reached the statistical equilibrium state (Boltzmann distribution). In that case, 
the system remains blocked in a quasi stationary state with a non standard distribution. We also study the 
relaxation of a test vortex in a steady bath of field vortices. The relaxation of the test vortex is described 
by a Fokker-Planck equation involving a diffusion term and a drift term. The diffusion coefficient, which 
is proportional to the density of field vortices and inversely proportional to the shear, usually decreases 
rapidly with the distance. The drift is proportional to the gradient of the density profile of the field vortices 
and is connected to the diffusion coefficient by a generalized Einstein relation. We study the evolution of 
the tail of the distribution function of the test vortex and show that it has a front structure. We also study 
how the temporal auto-correlation function of the position of the test vortex decreases with time and find 
that it usually exhibits an algebraic behavior with an exponent that we compute analytically. We mention 
analogies with other systems with long-range interactions. 

PACS. 05.20.-y Classical statistical mechanics - 05. 45. -a Nonlinear dynamics and nonlinear dynamical 
systems - 05.20.Dd Kinetic theory - 47.10.-g General theory in fluid dynamics - 47.32.C- Vortex dynamics 

1 Introduction of the Hamiltonian Mean Field (HMF) model, a toy model 

of systems with long-range interactions which possesses 

Systems with long-range interactions have been the ob- a lot of i nteresting properties and which can be studied 

ject of considerable interest m recent years (see different in great detaiL Qn the other handj Chavanis and 

contributions m the book [I]). Their dynamics is very rich Benedetti et al. [H] have worked out the kinetic theory of 

and presents many interesting features 03]. Therefore, a 2D C oulombian plasma and Valageas [T7] has built up a 

the construction of a kinetic theory adapted to such sys- kinetic theory for a m gravitat i 0Iia i system m a cosmolog- 

tems is a challenging problem with a lot of potential ap- icaJ context . Finally, Chavanis [T5] has obtained a general 

plications. Several kinetic theories have been developed kinctic equationj written in angle-action variables, which 

m the past. The first kinetic theory was constructed by is exp ected to describes the dynamical evolution of a one- 
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Indeed, it is not possible to assume that the system is 
spatially homogeneous as done in the kinetic theory of 
Coulombian plasmas (invoking Debye shielding [19j). in 
the kinetic theory of stellar systems El (invoking a local 
approximation [2D]), and for the HMF model (above the 
critical energy E c where a stable homogeneous phase ex- 
ists [10-12,14]). Furthermore, the point vortex gas is pe- 
culiar because point vortices have no inertia so that the 
collision term in the kinetic equation directly acts in po- 
sition space. Therefore, the role of position plays the role 
of velocity in usual kinetic theories. This makes the case 
of point vortices intermediate between kinetic theories of 
spatially homogeneous systems for which the evolution oc- 
curs only in velocity space [4-12,14-16] and kinetic theories 
of spatially inhomogeneous systems for which the evolu- 
tion occurs in velocity and position space |13j , or in energy 
space if we average over the orbits }17I18| . Because of its 
connections with other systems with long-range interac- 
tions, the kinetic theory of point vortices has a broader 
interest than simply fluid mechanics. 

A kinetic theory of point vortices in a shear flow has 
been first developed by Nazarenko & Zakharov [21] . They 
considered a multi-components system and assumed that 
the interaction between vortices is shielded due to geo- 
physical effects like rotation. They studied in detail the 
close collisions between two-screened particles (Stuart vor- 
tices) which are moved to each other by the collective 
shear flow and developed a kinetic theory a la Boltzmann. 
They showed that the mean vorticity profile does not 
change in time and that, due to collisions, the most in- 
tensive vortices are concentrated in the regions of large 
total vorticity while less intensive vortices are in small 
vorticity regions. For a single species system, the collision 
integral cancels out identically. 

A kinetic theory of point vortices with equal circula- 
tion has been developed more recently by Chavanis [22- 
24,2] (see also related works in Chavanis & Sire [25126] ') by 
analogy with the Brownian theory of Chandrasekhar [27] 
in stellar dynamics. In a first paper |22j . he considered the 
relaxation of a test vortex in a "sea" of field vortices and 
derived a Fokker-Planck equation where the evolution of 
the distribution function of the test vortex is due to the 
competition between a diffusion term and a drift term. 
The diffusion arises from the fluctuations of the velocity 
created by the field vortices and the drift term is due to the 
response of the field vortices to the perturbation caused by 
the test vortex, as in a polarization process. Chavanis |22j 
calculated the expression of the diffusion coefficient from 
the Kubo formula and the expression of the drift term 
from a linear response theory. This is similar to the calcu- 
lations of Kandrup |28j in stellar dynamics to determine 
the diffusion coefficient and the friction force experienced 
by a test star in a cluster of field stars. It is found that 



1 Of course, stellar systems are spatially inhomogeneous but 
the collision term is calculated as if the system were spatially 
homogeneous (this results in a logarithmic divergence of the 
diffusion coefficient at large scales). Then, spatial inhomogene- 
ity is taken into account in the inertial (Vlasov) term through 
the mean field gravitational potential. 



the diffusion coefficient of point vortices is proportional 
to the local vorticity created by the field vortices and in- 
versely proportional to the shear. On the other hand, the 
drift velocity is proportional to the local vorticity gradient 
and inversely proportional to the shear. Assuming that the 
vorticity profile of the field vortices is positive, axisymmet- 
ric and decreasing with the distance, the expression of the 
drift velocity shows that a test vortex with positive circu- 
lation climbs the vorticity gradient and that a test vortex 
with negative circulation descends the vorticity gradient. 
When the field vortices have the Boltzmann distribution 
of statistical equilibrium (thermal bath), the diffusion co- 
efficient and the drift coefficient (mobility) are related to 
each other by a generalization of the Einstein relation in- 
volving a negative temperature. In a second paper, Cha- 
vanis [23] developed a more complete kinetic theory of 
point vortices by using the projection operator formal- 
ism of Willis & Picard [29] . From this general formalism, 
he obtained a kinetic equation describing the evolution of 
the vortex system "as a whole". This is the counterpart 
of the generalized Landau equation obtained by Kandrup 
[50] in stellar dynamics. This equation conserves the en- 
ergy and monotonically increases the Boltzmann entropy 
(ii-theorem). The computed collision integral is of order 
0(1 /N) in a proper thermodynamic limit N — > +oo where 
the domain area V is fixed and the circulation of the point 
vortices scales like 7 ~ (so that the total circulation 
r = N~f is of order unity). This collision integral takes 
into account the influence of two-body correlations. For 
N — > +00, the correlations are negligible and we recover 
the 2D Euler equation which describes a "collisionless" 
evolution. This is the counterpart of the Vlasov equation 
in plasma physics. At order 0(1/N), the kinetic theory 
shows that the "collisional" evolution is due to a condi- 
tion of resonance between the trajectories of point vortices 
that can be far away. Therefore, this approach takes into 
account distant collisions between vortices while the ap- 
proach developed by Nazarenko & Zakharov [21] focuses 
on close collisions. This is why the collision integral de- 
rived in [23] can be non-zero for a single species system of 
vortices while it cancels out identically in [5T] . The above- 
mentioned kinetic theory developed at the order 0(1/N), 
describes the evolution of the system on a timescale of 
order Ntu, were to is the dynamical time. Furthermore, 
it assumes that the point vortices are transported by the 
collective shear flow (meanfield velocity) rather than, say, 
triple collisions. It is valid therefore when the shear is suf- 
ficiently strong in the system. If we implement a thermal 
bath approximation to describe the relaxation of a test 
vortex in a fixed distribution of field vortices, we recover 
the Fokker-Planck equation involving the terms of diffu- 
sion and drift obtained in [55]. In that case, the relaxation 
time scales like (N/ lnN)tr>. 

Similar problems have been studied independently by 
Dubin and collaborators [31-34] in the context of non- 
neutral plasmas under a strong magnetic field, a system 
isomorphic to the point vortex system. Dubin & O'Neil 
[5T] developed a kinetic theory of these systems, starting 
from the Klimontovich equation [35j and using methods 
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of plasma physics. Their kinetic theory is able to take 
into account collective effects that are ignored in the ap- 
proach of Chavanis [23| . However, the kinetic equation 
obtained by Chavanis [23] captures the main properties 
of the dynamics and can be solved numerically more eas- 
ily, an advantage that will be used in the present paper. 
In recent works, Dubin and collaborators [32-34] studied 
the dynamics of a test vortex in a background shear and 
provided very nice numerical simulations and laboratory 
experiments of this process. On a theoretical point of view, 
they derived the expressions of the diffusion coefficient and 
of the drift velocity of the test vortex. These expressions 
are consistent with those obtained previously by Chava- 
nis [22 23 . They also addressed the form of the cut-off 
to the logarithmic divergence that occurs in these quan- 
tities. However, they introduced these terms of diffusion 
and drift as independent effects Q and did not derive the 
Einstein relation connecting the diffusion to the drift nor 
the Fokker-Planck equation governing the evolution of a 
test vortex in a fixed distribution of field vortices. 

The object of the paper is to further develop the kinetic 
theory of point vortices and numerically solve the corre- 
sponding kinetic equations in order to illustrate the basic 
results of the theory. The kinetic theory is important to de- 
scribe the relaxation of the system towards the Boltzmann 
distribution predicted by equilibrium statistical mechan- 
ics for t — * +oo [36-41]. It is also necessary to determine 
the timescale of the collisional relaxation and to prove 
whether or not the system will truly relax towards Boltz- 
mann statistical equilibrium. Indeed, the convergence to- 
wards Boltzmann statistical equilibrium, which is based 
on an assumption of crgodicity, is not firmly established 
for complex systems with long-range interactions such as 
point vortices. A first reason is that such systems exhibit 
non-markovian effects and spatial derealization so that 
the monotonic increase of the Boltzmann entropy is diffi- 
cult to prove and could even be wrong in a strict sense |23j . 
It is only when additional approximations arc included 
(markovian approximation, neglect of three-body correla- 
tions,..) that the Boltzmann H-theorem is recover ed|. On 
the other hand, the evolution of the system is due to res- 
onances between different orbits and it may happen that 

2 The reason is that Schecter & Dubin [32] considered the 
drift of a test vortex in a shear created by a smooth vortic- 
ity field (without fluctuation) that is solution of the 2D Euler 
equation while Chavanis [22123] considered the drift of a test 
vortex in a shear produced by a discret collection of N point 
vortices. The fluctuations, due to finite N effects, give rise to 
the diffusion process and the inhomogeneity of the averaged 
vorticity profile gives rise to the drift. Since there are no dis- 
crete effects in the situation considered by Schecter & Dubin 
|32j . the test vortex experiences only a drift due to the inho- 
mogeneity of the vorticity background. Although the systems 
are different, the expressions of the drift obtained by Chavanis 
22 23] from the Liouville equation and by Schecter & Dubin 
|32| from the 2D Euler equation are the same. 

3 These approximations are expected to be valid in a proper 
thermodynamic limit TV — > +oo with fixed r\ = (3N^ 2 and 
e = E/N 2 j 2 . Therefore, the H-theorem holds in that limit 
N -> +oo. 



the evolution stops before the system has reached statisti- 
cal equilibrium because there is no resonance anymore. In 
that case, the system is blocked in a quasi stationary state 
(QSS) which can persist for a very long time. One object 
of the paper is to discuss these issues and perform numer- 
ical simulations in order to illustrate the particularity of 
the point vortex dynamics. 

The paper is organized as follows. In Seed] we briefly 
discuss the statistical mechanics of the point vortex gas at 
equilibrium and extend the theory to the case where the 
point vortices have different circulations. The predictions 
of equilibrium thermodynamics will be compared with the 
results of the kinetic theory throughout the paper. In Sec. 
[3] we develop the kinetic theory of point vortices and de- 
scribe the evolution of the system "as a whole" . In Sec. 13. H 
we generalize the kinetic theory of [23j , based on the pro- 
jection operator formalism, to a multi-species point vortex 
gas. Considering an axisymmctric evolution, we obtain an 
explicit kinetic equation for the evolution of the system at 
the order 0(1/N). We provide a simpler derivation of this 
kinetic equation than the one given in [23] . We also stress 
some technical difficulties associated with the kinetic the- 
ory and propose some solutions to circumvent these prob- 
lems. In Sec. I3.3[ wc show that the derived kinetic equa- 
tion conserves the vortex number of each species, the total 
energy and the total angular momentum. Furthermore, it 
increases the Boltzmann entropy (iJ-theorem) . For a sin- 
gle species system, the evolution of the point vortex gas 
on a timescale Nto is due to a condition of resonance 
Q(r,t) = f2(r',t) between different orbits r/r' that can 
be satisfied only when the profile of angular velocity is 
non-monotonic. As a result, the Boltzmann distribution 
corresponding to statistical equilibrium is not the only sta- 
tionary solution of the kinetic equation. For example, any 
distribution with a monotonic profile of angular velocity 
is a steady state solution. Therefore, the dynamical evo- 
lution stops when the profile of angular velocity becomes 
monotonic even if the system has not reached the Boltz- 
mann distribution (further evolution may occur on longer 
timescale due to terms of order 1/N 2 , 1/N 3 .... in the ki- 
netic theory corresponding to the influence of three-body, 
or higher, correlations). This kinetic blocking, described 
in Sec. 13.21 is illustrated numerically in Sec. 13.41 In Sec. [4] 
we consider the relaxation of a test vortex in a bath of field 
vortices. In Sec. I4.H we show that the stochastic process is 
described by a Fokker-Planck equation involving a term of 
diffusion and a term of drift. This Fokker-Planck equation 
can be derived from the projection operator formalism. 
Wc extend the approach of [23] by considering a test par- 
ticle with a circulation Pq that can be different from the 
circulation 7 of the field particles. In Sec. 14. 2\ we discuss 
specifically the terms of diffusion and drift. For a thermal 
bath of field vortices, we show that they are related to 
each other by a generalized Einstein relation £ = Df3To. 
Furthermore, the diffusion coefficient D(r) depends on the 
position and decreases rapidly with the distance, leading 
to anomalous dynamical behaviors. In Sec. 14. 3[ we derive 
the first and second moments of the radial position incre- 
ment from the Fokker-Planck equation. In Appendix [C] 
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we directly compute these moments from the Hamiltonian 
equations of motion. We can therefore justify the Fokker- 
Planck equation by this alternative approach instead of 
advocating the projection operator formalism. In Sec. 14.41 
we precisely discuss the scaling of the relaxation time with 
N. We distinguish the collisional timescale ~ Ntn associ- 
ated with the evolution of the system "as a whole" from 
the collisional timescale (N/ \nN)tjj associated with the 
evolution of a test vortex in a bath. In Sec. [5l we con- 
sider different features of the Fokker-Planck equation. In 
Sec. 15.11 we study the evolution of the tail of the distri- 
bution function of the test vortex and show that it has a 
front structure. We characterize the displacement of the 
front and its shape. In Sec. 15.21 we consider the tempo- 
ral auto-correlation function of the position of the test 
vortex and show that, in cases of physical interest, it de- 
creases algebraically with time. The exponent of algebraic 
decay is computed analytically and compared with direct 
numerical simulations. In Sees. [5J and [JJ we consider ex- 
plicit examples of Fokker-Planck equations corresponding 
to typical bath distributions of the field vortices. In order 
to be complete, we describe various examples of flow pro- 
files. However, we present numerical simulations only for 
the gaussian vortex. 

2 The point vortex gas at equilibrium 

2.1 The statistical equilibrium state 

Basically, the dynamical evolution of a system of point 
vortices in two dimensions is described by the Hamiltonian 
equations of motion [42] : 

dxj_ _ dH_ dyi_ _ dH_ 

ll dt ~ % ' ll dt ~ dx t ' u 

H = --^7i7jHri-r,-|, (2) 

»<; 

where 7, is the circulation of point vortex i. A particu- 
larity of this Hamiltonian system, first noticed by Kirch- 
hoff 43J, is that the coordinates (x,y) of the point vor- 
tices are canonically conjugate. The point vortex system 
conserves several quantities. The number N a of vortices 
of species a, or equivalently the total circulation of each 
species r a = N a "/ a , is conserved. The system also con- 
serves the energy E = H . There are additional conserved 
quantities depending on the geometry of the domain. The 
angular momentum L = ^ jirf is conserved in an infinite 
domain or in a disk and the linear impulse P = 7,rj 
is conserved in an infinite domain or in a channel. In an 
unbounded domain, we must consider vortices of the same 
circulation otherwise they would form pairs (dipoles) and 
escape to infinity, so there is no equilibrium. In a bounded 
domain, the Hamiltonian has to be modified so as to 
take into account the contribution of vortex images. Note 
that the Hamiltonian @ does not involve a kinetic term 
Y]j mjVf /2 usually present for material particles. This is 



because point vortices have no inertia. Therefore, a point 
vortex produces a velocity, not an acceleration (or a force), 
contrary to other systems of particles in interaction like 
electric charges in a plasma [5] or stars in a galaxy [5J. 
In a sense, this interaction is related to the conception of 
motion according to Descartes, in contrast to Newton. 

Assuming ergodicity, the point vortex gas is expected 
to achieve a statistical equilibrium state for t — > +00. 
The statistical mechanics of point vortices is very pecu- 
liar and was first discussed by Onsager [36] . He consid- 
ered box-confined configurations and showed that for suf- 
ficiently large energies, point vortices of the same sign 
tend to "attract" each other and group themselves to 
form "clusters" or "supervortices" similar to the large- 
scale vortices observed in the atmosphere of giant plan- 
ets. When all the vortices have the same sign, this leads 
to monopolcs (cyclones or anticyclones) and when vor- 
tices have positive and negative signs, this leads to dipoles 
(a pair of cyclone/anticlone) or even tripoles. These or- 
ganized states are characterized by negative tempera- 
tures. This is due to the fact that the structure function 
<P(E) = J H<E dx\dy\...dx]s[dyisi is finite for E — > +00 
since the phase space coincides with the configuration 
space: ^(+00) = J dxidyi-.-dx^dyN = V N where V is the 
area of the system. The existence of negative temperatures 
for point vortices should not cause surprise. For material 
particles, the temperature is a measure of the velocity dis- 
persion and it must be positive. Indeed, it appears in the 
Maxwell distribution e _/3u I 2 (or in the partition function) 
which can be normalized only for (3 > 0. However, since 
point vortices have no inertia, there is no such term in the 
equilibrium distribution of point vortices (see below) and 
the temperature can be negative. 

To obtain more quantitative results, Joyce & Mont- 
gomery [37j and Lundgrcn & Pointin [38] considered a 
mean field approximation. This is valid in a proper ther- 
modynamic limit N — > +00 in such a way that the nor- 
malized temperature 77 = /3N~f 2 and the normalized en- 
ergy e = E/N 2 j 2 are fixed (for a single species system). 
On physical grounds, it is reasonable to consider that the 
area of the domain and the total circulation of the vortices 
are of order unity. Then, by rcscaling the parameters ap- 
propriately, the proper thermodynamic limit corresponds 
to N -> +00 with 7 ~ 1/N, V ~ 1, E ~ 1 and /3 ~ N 
[2J. Since the coupling constant f3-f 2 — rj/N ~ \JN goes 
to zero for N — * +00, we are considering a weak long- 
range potential of interaction. In that limit, we can ne- 
glect the correlations between vortices and the statistical 
equilibrium distribution of point vortices is given by the 
Boltzmann distribution 

P eq (r) = Ae-WU, (3) 

where ip is the stream function. Using the Poisson equation 
Aip = —oj with the smooth vorticity w(r) = NjP(r), 
the statistical equilibrium state is obtained by solving the 
Boltzmann-Poisson equation 

- Alp = NjAe-^, (4) 
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and substituting the resulting stream function in Eq. ([3]). 
Joyce & Montgomery |37j obtained the equilibrium distri- 
bution j3]) by maximizing the entropy of the point vortex 
gas at fixed circulation and energy. Lundgren & Pointin 
[55] derived Eq. (fj| from the equilibrium BBGKY hierar- 
chy of equations in the limit N — > +oo. A rigorous deriva- 
tion of the mean field equations was provided by [39-41]. 

On the other hand, at sufficiently large negative en- 
ergies, the temperature is positive. In that case, like-sign 
vortices tend to "rcpell" each other. When all the vortices 
of the system have the same sign, they accumulate on the 
boundary of the domain. This regime can again be de- 
scribed by the mean field theory of [37-41]. Alternatively, 
when the system is neutral, opposite-sign vortices tend to 
"attract" each other resulting in a spatially homogeneous 
distribution with strong correlations between vortices. In 
that case, the point vortex gas is similar to a Coulombian 
plasma. This is the situation considered by Frohlich & 
Ruelle HI]. As discussed by Ruelle [35], we may expect a 
phase transition (related to the Kosterlitz-Thouless tran- 
sition) as a function of the temperature. At large (positive) 
temperatures, the system is in a "conducting phase" (in 
the plasma analogy) with free vortices that can screen ex- 
ternal "charges" . In that case, there can be Debye shield- 
ing like for a Coulombian plasma. At low (positive) tem- 
peratures, opposite-sign vortices tend to form pairs (+, — ) 
and the gas is in a "dielectric phase" where all charges are 
bound forming dipolar pairs. These "dipoles" are similar 
to "atoms" (+e, — e) in plasma physics. In that case, there 
is no screening. 

Summarizing, there are two very different regimes in 
the point vortex gas that correspond to different thermo- 
dynamic limits. In the mcanficld theory of [37-41], the 
thermodynamic limit corresponds to N — > +oo in such a 
way that the size V of the system is fixed and the inter- 
action between vortices is weak since the "charge" of the 
vortices tends to zero as 7 ~ 1/N. In that case, the sys- 
tem is spatially inhomogeneous, the correlations between 
vortices are negligible for N — ► +00, there is no shielding 
and the temperature can be negative. The physics of the 
problem is controlled by the one-body distribution func- 
tion P(r). Alternatively, Frohlich & Ruelle 44J consider 
a neutral system of 2N positive and negative point vor- 
tices with circulation ±7 and take the usual thermody- 
namic limit N, V -> +00 with N/V and E/N fixed (with 
7 ~ 1). In that case, the system is spatially homogeneous 
in average, the temperature is positive and the physics 
of the problem is controlled by the two-body correlation 
function. In the screening phase (at high positive tem- 
peratures), there is Debye shielding and the correlation 
function tends to zero exponentially rapidly. 

2.2 Generalization to a collection of circulations 

In this paper, we shall restrict ourselves to the situation 
described in [37-41]. In this section, we generalize the max- 
imum entropy method of Joyce & Montgomery [37] in 
order to determine the statistical equilibrium state of a 
system of point vortices with different circulations. This 



will be useful to interpret the results of Sec. 13.21 Follow- 
ing the Boltzmann procedure, we divide the domain into 
a very large number of microcells with size h 2 (ultimately 
h — ► 0). A microcell can be occupied by an arbitrary num- 
ber of point vortices. We shall now group these microcells 
into macrocells each of which contains many microcells 
but remains nevertheless small compared to the extension 
of the whole system. We call v the number of microcells 
in a macrocell. A macrostate is specified by the number 
{n-ia} of point vortices of circulation j a in the macrocell 
i (irrespective of their position in the cell) while a mi- 
crostate is specified by their precise position in the cell. 
Using a combinatorial analysis, the number of microstates 
corresponding to the macrostate {rii a } is 

W({n la })=l[N a \l[—. (5) 

a 1 

We introduce n a (r) the smooth density of point vor- 
tices of species a in r. The vorticity of species a is then 
w a (r) = 7a^ a ( r )- If w e define the Boltzmann entropy by 
S = InW, use the Stirling formula for n ia » 1 and take 
the continuum limit, we get 

S=-Y /^In^dr. (6) 

7a 7a 

The Boltzmann entropy ([6]) measures the number of mi- 
crostates corresponding to the macrostate specified by 
{w a (r)}. At statistical equilibrium, the system is expected 
to be in the most probable macrostate, i.e. the one that is 
the most represented at the microscopic level. Assuming 
that all the microstates are equiprobable, the equilibrium 
distribution is obtained by maximizing the Boltzmann en- 
tropy ([6]) while conserving the circulation of each species 
r a = J u) a dr and the mean field energy E = | J ujtpdr. 
The vorticity and the strcamfunction are related to each 
other by the Poisson equation 

-Ai) = u = Y,Ua- (7) 

a 

Finally, in an infinite domain or in a disk, we must also 
conserve the angular momentum L = j tor 2 dr. Introduc- 
ing Lagrange multipliers and writing the first order vari- 
ations as 

a a 5r a - f35E - X -iM L 8L = 0, (8) 

a 

we find that the most probable state is 

u a (r)=A a e-^'V, (9) 

where 

V = i> + \n L r 2 , (10) 

is the relative stream function. This describes a flow that 
is steady in a frame rotating with angular velocity fl^ (see 
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Appendix Therefore, the equilibrium state is obtained 
by solving the multi-species Boltzmann-Poisson equation 



-0Ja.il>' (r) 



(11) 



and substituting the resulting stream function back into 
Eq. ([9|). The Lagrange multipliers can then be related to 
the integral constraints. We note that the vorticity profiles 
of different species are related to each other by 



l/7a 



l/7b 



hence 



w„(r) = C of K(r)|W76 ; 



(12) 



(13) 



where C a b is independent on the position. Assuming that 
7a > 0, 7h > and that w (r) decreases with the distance 
(which corresponds to equilibrium states with f3 < 0), 
this relation indicates that intense vortices (7 a > 7b) are 
more concentrated at the center, on average, than weaker 
vortices. On the other hand, Eq. (fT3|) shows that opposite 
sign vortices tend to separate (at negative temperatures 
where the mean field theory applies). More generally, Eq. 
(fT3")) characterizes the seggregation between vortices with 
different circulations. 



3 Evolution of the system as a whole 

The equilibrium statistical mechanics tells nothing con- 
cerning the timescale of the relaxation towards statisti- 
cal equilibrium. It is furthermore not obvious whether the 
system of point vortices will truly relax towards statistical 
equilibrium or not. Indeed, the Boltzmann distribution Q 
is based on a hypothesis of ergodicity and on the assump- 
tion that all the accessible microstates are equiprobablc. 
This is essentially a postulate. In order to precisely an- 
swer these questions we need to develop a kinetic theory 
of point vortices. 



3.1 Kinetic theory of point vortices 

There are different methods to obtain a kinetic equation 
for a system of point vortices. One approach is to start 
from the Klimontovich equation and develop a quasi-linear 
theory as in plasma physics |31j . Another possibility is to 
start from the Fokker-Planck equation and calculate the 
first (Ar) and second moments ((Ar) 2 ) of the increments 
of position of a point vortex due to the interaction with the 
other vortices. This approach is developed in Appendix [Cl 
A third possibility is to start from the Liouville equation 
and use the projection operator formalism [23j . An inter- 
est of this formalism is to yield a general non-Markovian 
equation that is valid for flows that are not necessarily 
axisymmctric. The other formalisms assume from the be- 
gining that the distribution of point vortices is axisymmet- 
ric and work in Fourier space (for the angular variables). 



By contrast, the projection operator formalism remains in 
physical space which enlightens the basic physics. Indeed, 
even if the formalism is abstract and complicated [29] , the 
final kinetic equation takes a rather nice form which bears 
a clear physical meaning^. The drawback of that approach 
is that it ignores collective effects. For axisymmetric flows, 
these collective effects have been taken into account in the 
approach of Dubin & O'Neil [3Tj . 

We shall here extend the kinetic theory of Chavanis 
[23], based on the projection operator formalism, to the 
case of a multi-components point vortex gas. We shall not 
repeat the intermediate steps of the calculations that can 
be found in [23]. Generalizing these calculations in order 
to include a distribution of circulations among the point 
vortices, we obtain a kinetic equation of the form 

^ + (V) • v. Q = A £ f dr J dnVd - 0) 



xQ(t,t- t 



lbV u (l -> 0) 



d 



+7^(0-1) A 



dr v 

u) a (r, t - T)uj b (r 1: t - t), 



(14) 



where Q(t,t — r) is the Green function associated with 
the averaged Liouville operator constructed with the mean 
field velocity (V)(r, t) = — z x In words, this means 
that we must perform the time integration by moving the 
point vortices between t and t — t with the mean field 
velocity. On the other hand, V M (1 -> 0,t) = V^(l -» 
0,i) — (V^r, t) is the fluctuating velocity by unit of cir- 
culation created by particle 1 located in ri on particle 
located in r (see Appendix [B]) . To obtain this closed 
kinetic equation, we have implicitly assumed that three- 
body correlations are negligible. This neglect is valid at 
order 1/N in the proper thermodynamic limit N — > +oo 
defined above. In the general kinetic equation (TH)) , we 
clearly see the terms of diffusion and drift (first and sec- 
ond terms in the r.h.s.) and their connection to a general- 
ized form of Kubo formula (the integral over time of the 
velocity auto-correlation function). These points will be 
further developed in the sequel. The ratio of the "colli- 
sion" term (right hand side) on the advective term (left 
hand side) is of order 1 /N in the thermodynamic limit (see 
Sec. l4.4p . Therefore, this kinetic equation is valid at order 
0(1/ N) in an expansion of the collision term in powers 
of 1/N. Thus, it describes the evolution of the system on 
a timescale of the order Ntp where tjj is the dynamical 
time. 

We must distinguish two regimes in the evolution of 
the point vortex system: 



4 Recently, it has been found that this kinetic equation could 
also be obtained from a BBGKY-like hierarchy [46]. This con- 
siderably simplifies the formalism. The collision integral cor- 
responds to the term of order 0(1/ N) in a systematic expan- 
sion of the equations of the hierarchy in powers of 1/N when 
N — > +oo with fixed domain area V and fixed total circulation 

r = Nj. 
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(i) Collisionless regime: For fixed t and N — > +00, the 
collision term on the r.h.s. of Eq. (fT4|) vanishes and this 
equation reduces to the 2D Euler equation [23] : 



dt 



(V) • Vu a = 0. 



(15) 



This is the counterpart of the Vlasov equation in plasma 
physics [19] or stellar dynamics [20j (it is valid for each 
species of particles). It can be directly derived from the 
Liouville equation by neglecting correlations between vor- 
tices, i.e. by assuming that the TV-body distribution func- 
tion is a product of N one-body distribution functions [53] . 
The Euler equation (fT?j|) describes the collisionless evolu- 
tion of the point vortex system on a timcscalc t <C Ntjj 
where try is the dynamical time. For N 3> 1, the valid- 
ity of the Vlasov regime can be very long in practice. 
Starting from an unstable initial condition, the 2D Euler- 
Poisson system can undergo a "violent relaxation" driven 
by purely mean field effects [2] . This form of collisionless 
relaxation leads to the rapid emergence of a quasi station- 
ary state (QSS) on the coarse-grained scale [47-50]. 

(ii) Collisional regime: The "collision" term (right 
hand side in Eq. (|14p) takes into account finite N effects. 
It can be viewed as the first order correction of the Eu- 
ler/ Vlasov equation in an iV _1 expansion (see, e.g., [2] p. 
260). Therefore, the collision term in the kinetic equation 
(|14p manifests itself on a timescale ~ Nto- In this paper, 
we shall be essentially interested by this "slow collisional 
evolution" . 

In the following, we shall assume that the vorticity pro- 
file is a stable stationary solution of the 2D Euler equation, 
so that it evolves under the sole effect of collisions (finite 
N effects) on a timescale Ntjj. Therefore, we shall not 
describe the process of violent relaxation [47-50] due to 
mean field effects that takes place during the collisionless 
regime on a timescale to (see the Conclusion for a dis- 
cussion of the different regimes of the dynamics) . We thus 
start from a stable stationary solution of the 2D Euler 
equation (possibly resulting from a phase of violent relax- 
ation) and discuss its collisional evolution. If we restrict 
ourselves to axisymmetric flows, the kinetic equation (fl4| 
becomes 



76 K(1 0, t - t) — + la V ri (0 1, t - t) — 
or or 1 



xuj a (r,t- r)ub(ri,t- t), 



(16) 



where V r = V ■ f is the radial component of the veloc- 
ity in the direction f = r/r of particle (we adopt the 
same convention for particle 1). Since the collision term 
is of order l/N, the smooth distribution uj a (r,t) evolves 
on a timescale of order Ntjj at least, which is much larger 
than the timescale on which the fluctuations have essential 
correlations. Thus, we can make a Markov approximation 
cj(r, t — t) ~ w(r,t) and extend the time integral to +00 
(see, however, the end of this section). In the same order 
of approximations, we can consider that, to leading order 



in l/N, the point vortices follow circular trajectories given 

by 



»(i-r) = r, 6(t - t) = 6 - Q{r,t)T, 



(17) 



where (r, 9) specify the position of a point vortex at time 
t and fl(r,t) = (V)g(r,t)/r is the angular velocity cor- 
responding to the mean field velocity (V)e(r,t). Within 
these approximations, the kinetic equation (|16[) becomes 



Jb v r (i^o,t-T)— + 7a y ri (o^M-T)— 

or or 1 



where (see Appendix [Bl 
V r (l 0,t) 



n sm( 



2ir rf + r 2 — 2rr\ cos( 



(18) 



(19) 



K(i-o,t-<r) 



r x sin(0 - 0i - AQt) 



2tt rf+r 2 - 2rr x cos(6» - 6 1 - AQt) ' 

(20) 



with Afl = n(r,t) - f2(n,t). Since V ri (0 -> 1) = 
— (r/r i)V r (l — > 0), we can rewrite the foregoing equation 
in the form 



a 1 a r+00 r 2-K r +co 



15 Id 
xV r (l ->■ 0,t-T)[ 7fc-T7- - la — 7T- )w a (r,t)u) b (r\,t). 
r or n or 1 



(21) 



The kinetic equation involves the function 



M 



dexV r (1^0,t)V r (1^0,t-T). (22) 



Using Eqs. l[T9]) - (|20"]) . it is shown in [23|2] (see also Ap- 
pendix [B| that 



M 







1 - 









(23) 



where r< (resp. r>) denotes the smallest (resp. largest) of 
r and r\. Therefore, the final form of the kinetic equation 
is 



duj a 1 d ^ 



Ar dr 



r'dr'5(Q - f2') In 



1 - 


fe) 2 ] 







X Ib^b 



J 1 dUJg 

r dr 



1 du' b 

laUa-Jg-j), 



(24) 



where u>' stands for uj(r' , t) and Q' stands for Q(r', t). This 
generalizes the kinetic equation obtained by Chavanis [23] 
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for a single species of point vortices. The angular velocity 
and the vorticity are expressed in terms of the orthoradial 
velocity by 

dtb Id 
(V), = uj = --(r(V) s ), (25) 



so that 



c = l|-(r 2 /?). 

r or 



(26) 



For future convenience, it will be useful to write the log- 
arithm as 



hi 



E 

n=l 



lfr< 

n \r> 



(27) 



We note that, for a multi-species system, a logarithmic 
divergence J2n=i 1 / n a PP ears m the kinetic equation l[2"l]) 
when r' = r and b ^= a. This problem will be discussed 
specifically in Sec. 14.11 We shall hcuristically regularize 
the divergence by introducing an upper cut-off A in the 
series ([2"T]) when necessary, writing J2n=i n( r </ r >) 2 ™- 

We also emphasize that the Dirac function arising in 
Eq. (|24[) does not make sense for values of r and r' such 
that Q(r) = Q{r') and (df2/dr)(r') = (in that case, 
the identity (|5"2"|) clearly breaks down). Mathematically 
speaking, the "function" x — ► S(x 2 ) has no sense even in 
the space of distributions. To overcome this problem, we 
can notice that this Dirac distribution comes out formally 
from Eq. (|2"2"|) where the time integration takes over all the 
interval [0, +oo). Therefore, one could expect that if the 
time averaging is taken on a finite interval [0, t) only, as in 
the initial integral ([TB]) . then the resulting kernel M(t) will 
be well defined as a smoothing approximation to Eq. (|23l) . 
The corresponding computation is given in Appendix B 
and leads to 



M(t) 



1 1 

47rr 2 AO 



arctan 



R 2 sm(tAC2) 



1 - R 2 cos(tAf2) 



(28) 



where R = r</r>. For t — > +oo, we recover formula ()23D . 
However, the regularized expression ([28]) does not suffer 
the problem discussed above (we may also wonder whether 
collective effects [21] that we have neglected can regularize 
or not the integral in the situation mentioned above). 



3.2 Condition of resonance and kinetic blocking 

For a single species of point vortices, the kinetic equation 
(J2U) becomes [23]: 



4r dr 



1 - 









, 1 duj 1 duj' 

x i uj _— - w — _ 

r or r or 



(29) 



This kinetic equation is the counterpart of the Lan- 
dau equation describing the dynamical evolution of a 3D 



plasma or a stellar system as a whole [30]. As is well- 
known, the Landau equation yields a logarithmic diver- 
gence J +O ° dk/k at small and large scales. The small scale 
divergence is regularized at the Landau length, corre- 
sponding to a deflection at 90° of the particles' orbits so 
that the linear trajectory approximation made by Landau 
is not valid anymore. The large scale divergence is reg- 
ularized, in plasma physics, by the Debye shielding (as 
shown by the Lenard-Balescu treatment of collective ef- 
fects) and, in stellar dynamics, by the finite size of the 
system (Jeans length). For the single species point vortex 
gas, we stress that, contrary to the Landau equation, there 
is no logarithmic divergence in the kinetic equation (1291) . 
Therefore, although there is no Debye shielding in the sin- 
gle species point vortex gas, the kinetic equation ([29]) is 
well-posed mathematically 0. As discussed in [3112312] . the 
"collisional" evolution of point vortices described by Eq. 
([29]) is due to a condition of resonance which can be satis- 
fied only if the profile of angular velocity is non-monotonic 
0. The current of vorticity in r is due to long-range colli- 
sions with point vortices in r' ^= r whose orbits satisfy the 
condition 0(r,t) = 0(r',t). The self- interaction at r = r' 
does not produce transport since the term in parenthesis 
vanishes identically. When the profile of angular velocity 
becomes monotonic, the evolution stops and the system 
becomes "frozen" in a QSS satisfying dui/dt = 0. This 
QSS usually differs from the statistical equilibrium state 
as will be shown numerically in Sec. 13.41 In that case, the 
relaxation towards statistical equilibrium (if it really hap- 
pens) takes place on a timescale larger than Nto which is 
not described by the present approach. If we want to de- 
scribe this regime, we need to take into account terms of 
order l/N 2 or smaller in the kinetic theory. They are as- 
sociated with three-body (or higher) correlation functions 

|23| . 

The situation is different if the system consists in a 
collection of point vortices with different circulations. As- 
suming that the profile of angular velocity is monotonic 
and using 5(i7 — Q') = S(r — r')/\^-(r)\, the kinetic equa- 
tion (EH) becomes 



duj a 

at 



i 

4r 



d_ 

dr 



h^-E 



i 



^r,t)\ 



dr 



dr 



JaUJa 



duJb\ 

dr y 

(30) 



where a cut-off A has been introduced in Eq. ([27]) so that 
In A ~ X)n=i We see that the diffusion current per 



5 The situation is different for a unidirectional flow where 
a logarithmic divergence occurs at large scales for the single 
species point vortex system (see Appendix E.2. of [23]). In 
that case, it must be regularized by invoking the finite extent 
of the system (like for self-gravitating systems) or geophysical 
effects like the finite Rossby radius of deformation that plays 
a role similar to the Debye length in plasma physics. 

6 We recall that, within our assumptions, this non-monotonic 
angular velocity profile must be stable with respect to the 2D 
Euler equation so as to avoid a "violent relaxation" process 
driven by mean field effects. It is indicated in [31] that such 
profiles can be realized experimentally [51] . 
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species does not vanish anymore when there arc at least 
two species of particles in the system Q. In that case, the 
kinetic equation (|30|) becomes very similar to the one ob- 
tained by Nazarenko & Zakharov [2T 0. The transport 
of point vortices of species a is caused by collisions with 
point vortices of other species at the same radial distance 
r = r' . However, using the anti-symmetry of the collision 
term, we can easily see that the global vorticity distribu- 
tion a/ = uj a does not change with time, i.e. 



(31) 



Finally, using the ii-theorem of Sec. 13.3.41 it is simple 
to show that the stationary solution of Eq. ([3"0|) (for all 
species) corresponds to 



du) a duj b 
or or 

for any a and b. This is equivalent to 

w.(r) =C ab \uj b (r)\^/^, 



(32) 



(33) 



where C a b is independent on r. This relation was previ- 
ously obtained in [21134] . It is similar to the relation Q13p 
derived for point vortices at statistical equilibrium. From 
Eq. (f32|) or (|33[) , we find that the vorticity of each species 
can be written 



w (r)=A,e- /57 " x(r) , 



(34) 



where A a and (3 are some constants and x( r ) is determined 
by the initial conditions. However, contrary to the case of 
statistical equilibrium © considered in Sec. 12.21 x(r) does 
not represent in general the stream function |34j . 



3.3 Conservation laws and H-theorem 

In this section, we show that the kinetic equation (|24|) re- 
spects the conservation laws of the point vortex dynamics 
and increases the Boltzmann entropy (iZ-theorem). We 
write the kinetic equation as 



dt 



ldJa 

r dr 



(35) 



where 



1 r + °° 
J a = / r'dr'S(0 - Q 1 ) In 

L J 







1 - 









)_dc4 

' r' dr' 



7 As shown by Dubin [31] , the kinetic theory is only valid for 
retrograde vortices such that (j a + *y b )E < 0. The sum in Eq. 
(|30l) must take into account this constraint. 

8 The prefactor in front of the parenthesis is, however, dif- 
ferent because we are not modelling the collisions exactly in 
the same way. We are considering collisions with large impact 
parameters while Nazarenko & Zakharov [21] consider close 
binary collisions; see the discussion in the Introduction. 



denotes the diffusion current. Due to the conservative 
form of Eq. ([35]) . it is clear that the total circulation 
r a = J uj a dr of each species of point vortices is conserved 
provided that J a vanishes at the fronticrc of the domain. 



3.3.1 Boltzmann distribution 

The Boltzmann distribution 

u a =A a e-^', (37) 

is a stationary solution of the kinetic equation (|24p. In- 
deed, using 

duj a „ dtp' ( 'dip 



(3 la uj a r(n -Q L ), (38) 



we find that the term in parenthesis in Eq. p6|) is equal 
to 



, 1 duj a 1 dui'j . 

lb^b~-^ laUa—.-Z-T = PlalbUaU b {Sl ~ 11 )• 



t 1 dr' 



(39) 



Therefore, the integrand in Eq. (|36p is proportional to 

5{Q- Q'){Q- Q') = 0. (40) 

Therefore, the current vanishes and the Boltzmann dis- 
tribution is a stationary solution of the kinetic equation 
([2"4")l . However, this is not the only stationary solution. For 
a single species system, any vorticity distribution with a 
monotonic profile of angular velocity is a stationary solu- 
tion of the kinetic equation ([2"9"]> since f?(r) 7^ for 
any r 7^ r' (and for r — r' the term in parenthesis van- 
ishes). For a multi-species system, the steady distributions 
of the kinetic equation (|24[) with a monotonic profile of an- 
gular velocity are given by Eq. (|34|) . They are in general 
different from the Boltzmann distribution. 

3.3.2 Conservation of energy 

The time variation of energy can be written 



E 



1 9u} j 



E 



ib - a 2irrdr 
at 



= 2nY J a -fdr ^ -2ttY J a Qr dr. (41) 

To get the first equality, we have used the Poisson equa- 
tion (J7J) and integrated by parts twice. To get the third 
equality, we have used Eq. ([35]) and integrated by parts. 
Inserting the current ([36]) in Eq. (|41|). we get 



E = — 



E 

a, b 



-foe 



rr'drdr'5(n- [2')fl\n 

I 1 duj a 







1 - 
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Interchanging the dummy variables a, b and r, r' , we ob- 
tain 



7T f + oc 
E=-^2 rr'drdr'5(Q - Q')Q' In 

a,b J° 

1 1 du a 



1 



1 duji 



<i^-^r-7^ v ^. (43) 

Taking the half-sum of these two expressions, we find that 
E = -^Yj / rr'drdr'8{n - Q')(Q - Q') 

a.b 

Xl T-(^jJHr^-^V^J- (44) 
Using Eq. (|4*U)) . we conclude that -E = 0. 

3.3.3 Conservation of angular momentum 

The time variation of angular momentum can be written 



L 



du) a 



Of 

" + 00 

/() 



2irrdr 



C 

/•-t-oo 

= 4tt V / J a rdr. (45) 
Inserting the current (f3"6"| in Eq. (j4"5)) . we get 

/■+00 

L = — 7r N / rr' drdr' S(il — f2') In 

: r jo 



a,b 



j 1 (9lJ a 

r <9r 



1 dul 

di 



(46) 



Interchanging the dummy variables a, b and r, r' , we ob- 
tain 



L= ^L 

i Jo 

a.b 



rr'drdr'8(fl - Q') In 



1-1^ 



2 1 



, 1 8uJa_ _ 

b r dr ^ a a r' dr' 



ldu' b 



(47) 



Taking the sum of the last two expressions, we conclude 
that L = 0. 



3.3.4 H-theorem 

The time variation of the entropy ([6]) can be written 

" + °° 1 diO a 



„ Jo 



-fa^a dr 



J a dr. 



(48) 



To get this expression, we have used Eq. (|35[) and inte- 
grated by parts. Inserting the current (|3"6"1) in Eq. (|4"8")h we 
get 



x In 



1 



a, 6 

'■> 



+°° i a > 

drdr'—r'^5(n-n') 

JaLUa Or 

,lduja__ 1 duj' b 
r or r or' 



(49) 



Interchanging the dummy variables a, 6 and r, r' , we ob- 
tain 

s = ^V [ + °° drdr > J- r ^A S (f2-n>) 



x In 



a, l> 

1^ 



, 1 3w a 



1 dw£ 



r' dr' 



Taking the half sum of these two expressions, we find that 



S 



x In 



1 

rr drdr 8(Q — fl ) 



1 - 



dr 



r' <9r' 



(51) 



Since this quantity is positive, the entropy cannot de- 
crease. An ff-theorem results: S > 0. However, the con- 
dition 5 = does not only select the Boltzmann distribu- 
tion. It is satisfied for any steady solution of Eq. (124]) . 

3.4 Numerical simulations 

We have performed numerical simulations of the kinetic 
equation (f!?9")) for a single species system of point vortices. 
To solve this equation, we have used the identity 



8{Q{r)-Q{r'))=Y, 



8(r' - r k ) 
\mrk)\ 



(52) 



where {rk} is the set of resonant points which satisfy the 
condition f2(r k ) = Q(r). Of course, this relation is valid 
only for resonant points r k such that 



dfl 
dr 



(r*) + 0. 



If a resonant point with Q!(r k ) = is encountered, then 
the kernel of the collision kinetic operator has to be mod- 
ified and one of the many different ways to regularize its 
expression is presented at the end of section 3.1 (it may 
also be necessary to reconsider the approximation (fTT)) 
when f2'(r k ) = 0). Assuming here that this situation never 
occurs, the kinetic equation (j2T))) becomes 



dco 
~dt 



1 



'^dr^\™(r k ,t)\ 









In 


1 - 











. . r k duj dui . . 

U)(r k ,t) — - uj—(r k ,t) 

r or or 



(53) 
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where r< (resp. r>) denotes the smallest (resp. largest) of 
r and r k . This can be written in the more compact form 



duj 
~dt 



ld_ 

r dr 



r[D(r,t) 



with 
D(r,t) 



.JLV — 



rk 



(r k ,t)\ 



In 



3a; 
<9r 



In 



ujr/(r, t) 



(54) 







I— 1 

1 









(55) 



'^(r*,*)l 



<9r 







I— 1 

1 









9a; 
dr 



(r k ,t). 



(56) 



We have considered an initial vorticity field iOo(r) as- 
sociated with a non-monotonic profile of angular veloc- 
ity fioif)- Specifically, we have taken ujo(r) = cxp[— (r — 
0.2) 2 /0.05] in a disk of radius R = 1 (see Figs. Q] and [2}. We 
numerically find that the system evolves until the profile 
of angular velocity becomes monotonic (see Fig. ^ . In the 
present situation, the vorticity profile also becomes mono- 
tonic (see Fig. [J}. When there is no resonance, the evolu- 
tion stops and the system remains blocked in a stationary 
state generically different from the Boltzmann distribu- 
tion 0. Therefore, the main effect of long-range collisions 
between point vortices is to make the profile of angular ve- 
locity monotonic. To our knowledge, the kinetic theory of 
point vortices |3f I23j is the first kinetic theory to exhibit 
such a behavior. In usual kinetic theories developed for 
ordinary gas [J], plasmas [5] and stellar systems [BJ, the 
system described by the Boltzmann equation (or by the 
Landau or Lenard-Balescu equation) always relaxes to- 
wards the Boltzmann distribution (the case of stellar sys- 
tems is peculiar because of the phenomena of evaporation 
and gravothermal catastrophe [55] but the convergence to 
the Boltzmann distribution is, however, the general ten- 
dency). Alternatively, for the homogeneous phase of the 
HMF model (or other one-dimensional systems) , the colli- 
sion term vanishes identically |lf I12I3| at the order f /N so 
that there is no evolution at all on a timescale ~ Nto (nu- 
merical simulations [53] show that the HMF model relaxes 
towards statistical equilibrium on a timescale N 1 ^to). 
For point vortices, the situation is intermediate between 
these two extremes. The system described by the kinetic 
equation (|2T)]) evolves until the profile of angular velocity 
becomes monotonic, then stops. For the inhomogencous 



9 This has been checked numerically by computing the term 
in parenthesis in Eq. (|29p for different couples of points r and 
r' . This term is not proportional to u(r)u(r')[i?(r) — 0(r')] 
as would be the case for a Boltzmann distribution according 
to Eq. ((39}. Note that the tail of the distribution does not 
evolve with time (since vortices in the tail are never in res- 
onance with vortices in the core) so that the tail is clearly 
non-Boltzmannian. The previous check shows that the final 
distribution of the core (that has evolved through resonances) 
is not Boltzmannian neither. 



0.75 




0.25 



Fig. 1. Evolution of the vorticity profile obtained by solving 
numerically the kinetic equation (|29|l . This kinetic equation is 
valid on a timescale Nto- On this timescale, the vortex gas 
does not reach statistical equilibrium but remains blocked in a 
QSS with a monotonic profile of angular velocity (see Fig. [2]| . 




Fig. 2. Evolution of the profile of angular velocity correspond- 
ing to the vorticity field represented in Fig. [JJ The evolution 
stops when the profile of angular velocity becomes monotonic 
so that there is no resonance anymore. Note, parenthetically, 
that the terminal point of the curve does not move. Indeed, 
according to Eq. (|125[) . the angular velocity at r = R is 
fl = r/(2irR 2 ) which is constant. 



phase of the HMF model (or other one-dimensional sys- 
tems), we expect a behavior similar to that observed for 
point vortices when we use angle-action variables, as dis- 
cussed in [18]. However, the problem is more difficult to 
investigate numerically because there is a richer variety of 
resonances. 

Let us conclude this section by some remarks that 
should be given further consideration in future works: 

(i) The diffusion current vanishes at r either because 
there is no point r' ^ r such that f2(r') = Q{r) (no res- 
onance) or because the term in parenthesis vanishes for 
any r' such that fl{r') = Q{r). Therefore, we have three 
possibilities: (i) the profile of angular velocity is mono- 
tonic everywhere so that the current vanishes because of 
the Dirac function (ii) the vorticity profile is the Boltz- 
mann distribution so that the current vanishes because 
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the term in parenthesis is zero at each point of resonance 
(if any) (iii) we are in a mixed situation where there is no 
resonance in certain regions and the term in parenthesis 
is zero at each point of resonance in other regions. 

(ii) As mentionned above, formula (|52[) is wrong when 
Si 1 {r^) = 0, i.e. when r is in resonance with a point lo- 
cated at an extremum of angular velocity 0. This means 
that initial conditions that have at least two extrema (in 
addition to r = 0) are not allowed by this model. Further- 
more, it is not proven that initial condition having at most 
one extremum will not generate profiles with two extrema 
in finite time (although this seems relatively unlikely). In 
other words, it is not clear whether the model will break 
down or not at finite time for some initial conditions. To 
overcome these problems, a smoothing strategy has been 
proposed at the end of section 3.1, and developed in Ap- 
pendix B. 

(iii) If the profile of angular velocity presents zones of 
resonances and zones of non-resonance, then the diffusion 
coefficient (|55[) is discontinuous at the separation and this 
can lead to numerical and physical instabilities. 

Therefore, the complete study of the kinetic equation 
(|29|) is complicated. We have presented here just one ex- 
ample of evolution. A more thorough (numerical and theo- 
retical) study of this equation would be certainly valuable. 
We also recall that the results of this section assume that 
the vorticity profile ujq(t) is stable with respect to the 
2D Euler equation. Now, a system with a non-monotonic 
profile of angular velocity can be unstable to diocotron 
modes [51]. This diocotron (or Kclvin-Hclmholtz) instabil- 
ity, which is a non-axisymmetric instability, can develop in 
the nonlinear regime and lead to a violent collisionlcss re- 
laxation [55-57]. The evolution is therefore very different 
from the one reported in Figs. Q] and [2] However, there 
also exists systems with non-monotonic profile of angu- 
lar velocity that are stable with respect to the 2D Euler 
equation |51|31] and that evolve under the sole effect of 
"collisions" like in Figs. [T] and [5] A better characterization 
of these regimes (and their selection) would be valuable. 



that we wish to describe. This problem is simpler than the 
previous one (evolution of the TV-body system as a whole) 
because it assumes that the distribution of the field par- 
ticles is fixed, while in the preceding case it was evolving 
self-consistently. This fixed distribution either represents 
the statistical equilibrium state (thermal bath) which does 
not change at all or a quasi-stationary distribution that 
evolves on a timescale that is much larger than the typical 
relaxation time of the test particle. We shall now consider 
this classical problem in the context of two-dimensional 
point vortices. 



4.1 The Fokker-Planck equation 

We consider the relaxation of a test vortex with circulation 
To in a "sea" of field vortices with circulation 7 (the field 
vortices play the role of the bath). In the TV — > +00 limit, 
the TV-body distribution of the field vortices can be written 



I'h 



(57) 



where Pbath(r) is the steady probability distribution of 
a single vortex of the bath. The vorticity profile of the 
bath is uj(r) = NjPb at h(r). Using the projection operator 
formalism, we find that the evolution of the density prob- 
ability P(r, t) of finding the test vortex in r at time t is 
governed by an equation of the form |23j : 



dP d /"* f 



xg(t,t-r) 
-rbv^o 



7V(1 -» 0,t 



0,i) 

)— 



d 
~dr[ 



dr v 

P(r,t-r)w(ri). (58) 



For an axisymmetric flow, we can repeat the same steps 
as in Sec. 13. II and we successively obtain 



4 Relaxation of a test vortex in a sea of field 
vortices 

Another classical problem in kinetic theory concerns the 
description of the relaxation of a test particle in a bath 
of field particles. If we inject a test particle in a steady 
bath of field particles, it will undergo a stochastic process 



10 In the numerical simulation reported in Figs. [T] and [5] the 
diffusive nature of the kinetic equation implies that ui'(0,t) — 
Q'(Q,t) — for t > (even if this is not the case initially). 
Therefore, at each time t > 0, there exists a point r > which 
is in resonance with the point = at which f2'(rk,t) = 0. 
However, the point rk = is special. Indeed, the numerator 
of D(r,t) is proportional to rfcln(l — (rk/r) 2 ) ~ r\ and the 
numerator of r?(r, t) is proportional to ^'(r^, t) ln(l — (r^/r) 2 ) ~ 
r\. These terms compensate the term f2'(rk,t) ~ — > in the 
denominator of Eqs. (|55|l - (|56|l . so that there is no singularity: 
D(r,t) ~ r){r,t) ~r|/r fc 



0. 



dP_ 

~dt 



1 d 



r dr 
xV r (l 0,t 



X) />27T 

dr / dO 

Jo Jo 

\( 1 8 r 
T ) 7-— - "0 



rr\driV r {l 
1 d 



d 



r 



7"i dr\ 



0,t) 

p{r, tMn), 



(59) 



and 



dP 
~dt 



^- = —4- / r'dr'5(f2 - St) In 
4r dr 1 



+00 



1 - 









,ldP p n ldw' 
r or r dr' 



(60) 



Equation ([50]) can be viewed as a Fokker-Planck equation 
(see below) describing the relaxation of the test vortex. 
It can be obtained directly from Eq. ([M]) by replacing Lo a 
by the distribution of the test particle P(r, t) and to>6 by 
the static distribution of the bath to(r). This procedure 
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transforms the integrodifferential equation (f24|) into a dif- 
ferential equation ((60|) . 

Let us first assume that the field vortices are at sta- 



tistical equilibrium so that w(r) = Ae 



represents 



the Boltzmann distribution (where ip' denotes the rela- 
tive stream function). Introducing 



1 dw 
r' dr 



(r') 



1 



" J — 

r' ' dr 



(61) 



in Eq. ([pTJ]) . using the S- function allowing to replace Q' by 
fi, and using Eq. (|6ip again with r instead of r', we find 
that the Fokker-Planck equation (|60|) can be rewritten in 
the form 



dP 
~dt 



ld_ 

r dr 



rD(r) 



dP 
dr 



(3r Q p 



dtp' 

dr 



with a diffusion coefficient 



D(r) 



7 
4r 2 



r'dr'5(f2- Q')ln 



o 



(62) 



;(/)• 
(63) 



These equations are valid even if the profile of angular ve- 
locity is non monotonic since the Boltzmann distribution 
is always a steady state of the kinetic equation (f2~4")l . This 
statistical equilibrium distribution does not evolve at all. 
We see that the diffusion coefficient in r is due to interac- 
tions with vortices whose orbits satisfy f2(r') = Q(r). This 
includes the local interaction with vortices at r' = r but 
also the interactions with far away vortices with r' =/= r. 
To our knowledge, this is the first kinetic theory where the 
diffusion coefficient exhibits such a spatial derealization. 

We now consider a bath with a monotonic profile of 
angular velocity that is not necessarily the Boltzmann dis- 
tribution of statistical equilibrium. Indeed, we have seen 
in the previous section that any Euler stable distribution 
with a monotonic profile of angular velocity is a station- 
ary solution of the kinetic equation (j29|) . Therefore, this 
profile docs not evolve on a timescale of order Nto on 
which the kinetic equation (|29|) is valid (but it may evolve 
on a longer timescale). We shall see that Ntjj is precisely 
the timescale controlling the relaxation of the test vortex, 
i.e. the time needed by the test vortex to acquire the dis- 
tribution of the bath. Therefore, we can consider that the 
distribution of the field vortices is frozen on this timescale. 
Thus, let us consider a bath that is not necessarily a ther- 
mal bath but that is a slowly evolving out-of-equilibrium 
distribution. If the profile of angular velocity of the bath 
is monotonic, then 



(64) 



In that case, the Fokker-Planck equation ((60|) can be writ- 
ten 



dP 1, ,13 
— = - \iiA- — 

at 4 r or 



I dP _, ^du) 

-yuj— — r P— 



\S(r)\\' dr 



dr 



, (65) 



where S(r) = ri?'(r) is the local shear created by the field 
vortices and wc have introduced the Coulomb factor 



\nA = Y -. 
^ n 

n=l 



(66) 



We note that the series diverges when n — > +oo. This di- 
vergent sum occurs because nearby vortices following un- 
perturbed orbits take a long time to separate. However, 
our theory breaks down at small separation because for 
separations less than d (say) we cannot assume that the 
motion of the particles is given by the unperturbed trajec- 
tory l|17[). In that case, one must consider the detail of the 
interaction between neighboring vortices (note that a sim- 
ilar treatment is required in 3D plasma physics and stel- 
lar dynamics to regularize the logarithmic divergence of 
the diffusion coefficient at small scales, i.e. at the Landau 
length [7]). Phenomenologically, the logarithmic diver- 
gence can be regularized by adding cutoffs so that In A = 
h\(r/d) |22|23j . A precise estimate of the lower cut-off d 
has been given by Dubin and collaborators [32-34]. They 
propose to take d = Max(<5, I) where I is the trapping dis- 
tance I = (2 7 2 /[47r(7 + lb)!!?!]) 1 / 2 and S is the diffusion- 
limited minimum separation S = (4D/|Z'|) 1 / 2 (where D 
is the diffusion coefficient given by Eq. (|69|) ). Orders of 
magnitude indicate that 5/1 ~ (-D/7) 1 / 2 ~ (In A) 1 / 2 and 
r/l ~ i?[( 7 + To^l/V] 1 / 2 ~ i?[( 7 + r )Nj/j 2 R^ ~ 
[(7 + r , o)/7] 1/2 iV 1/2 . Therefore, the logarithmic factor 
scales with N as 



hi A 



-kiN, 
2 



(67) 



in agreement with the rough estimate given in [2"5] . 

The Fokker-Planck equation (|65"|) can be rewritten in 
the form 



dP 
~dt 



ld_ 

r dr 



. dP r dlnU 

#h -p— 

or 7 dr 



with a diffusion coefficient 

7 



D(r) 



In A- — — — 

4 \E(r)\ 



(68) 



(69) 



This is a drift-diffusion equation describing the evolution 
of the test vortex in an "effective potential" 



ET eff (r) = -(r /7)lnKr)| ) 



(70) 



produced by the field vortices [23] . The diffusion coefficient 
is proportional to the density of field vortices to(r) and 
inversely proportional to the local shear |i7(r)| created 
by the background vorticity distribution. The equilibrium 
distribution of the test vortex is 



P e (r)cx \uj{r)\ r °'\ 



(71) 



up to a normalization factor. When the test vortex has the 
same circulation as the field vortices (Jo = 7), the test 
vortex ultimately acquires the distribution of the bath: 
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P e (r) oc |w(r)|. However, when the test vortex has a cir- 
culation different from that of the bath, the distribution of 
the test vortex differs from the distribution of the bath by 
a power r Q /-f. If the field vortices are at statistical equi- 
librium (thermal bath), their distribution is determined 
by the Boltzmann-Poisson equation 



u(r) = -Aip = Ae 



(72) 



Assuming that the profile of angular velocity is monotonic, 
we can use the Fokker-Planck equation ([68]) with In \ u>\ = 
-/?7V>' + In \A\ to obtain E3 



dP 
~dt 



ld_ 

r dr 



(73) 



with a diffusion coefficient given by Eq. ([69]) . The Fokker- 
Planck equation (|73|) describing the relaxation of a point 
vortex in a sea of field vortices is the counterpart of the 
Kramers equation derived by Chandrasekhar [27] to de- 
scribe the relaxation of a test star in a cluster. We em- 
phasize that the diffusion coefficient in Eq. (|73p depends 
on the position r of the test vortex [22]. Similarly, in the 
Kramers-Chandrasekhar equation, the diffusion coefficient 
of the test star depends on its velocity v [27] . 



4.2 Diffusion coefficient and drift term 

We can write the Fokker-Planck equation ([59]) in a form 
that explicitly isolates the terms of diffusion and drift: 



dP 

"at 



r dr 



dP 
D—-P V 
or 



(74) 



The diffusion coefficient is given by 

r+oo r 

D = lJ dT J drxV r {l -► 0,i)K(l 0. 



t - r)u(ri). 



It can be written in the form of a Kubo formula 



D = 



(V r (t)V r (t - T)) dT, 



(75) 



(76) 



representing the time integral of the velocity auto- 
correlation function [23j . The drift term is given by 



v = r 



dr / driK(l -> 0, t)V r (l -»• 0,t - r) 



r du>i 
r\ dr\ 
(77) 



11 We note the analogy of Eq. (|73p with the ordinary 
Smoluchowski equation describing the sedimentation of col- 
loidal suspensions in a gravitational field. Usually, the Smolu- 
chowski equation is obtained from the Fokker-Planck equation 
(Kramers equation) in a strong friction limit where the inertia 
of the particles can be neglected. In the present context, since 
the point vortices have no inertia, the Fokker-Planck equation 
describing the relaxation of a test particle directly has the form 
of a Smoluchowski equation in physical space. 



It can be seen as a sort of generalized Kubo formula in- 
volving the gradient of the density of field vortices instead 
of their density itself. The physical origin of the drift can 
be understood by developing a linear response theory |22j . 
It arises as the response of the field vortices to the per- 
turbation caused by the test vortex, as in a polarization 
process. 

According to Eq. ([6TJ]) , the diffusion coefficient and the 
drift term can be written more explicitly as 



D 



7 

4r 2 



r'dr'5{Q - Q') In 







1 - 









p r+oo 

— / dr'8(Q-Q')hx 
4r Jo 



r> 



duj' 
dr 7 ' 



(78) 



(79) 



When the profile of angular velocity is monotonic, using 
Eq. ([64]) , we find that the diffusion coefficient is given by 



D(r) = Jln^ W (r). 



(80) 



This expression of the diffusion coefficient, with the shear 
reduction, was first derived in Chavanis [22|23j . An equiv- 
alent expression has been obtained by Dubin & Jin [33] . 
On the other hand, using Eq. (|64[) . we find that the drift 
term is given by 



'/ 



-In^^L 

4 \S(r) 



du> 



dr 



(r). 



(81) 



The drift is proportional to the local vorticity gradient 
and inversely proportional to the shear. Comparing with 
Eq. ([50]) . we find that the drift velocity is connected to 
the diffusion coefficient by the relation 



n 



7 dr 



(82) 



This relation is valid for an arbitrary distribution of field 
vortices provided that the profile of angular momentum 
is monotonic. This expression was given in Chavanis [23] 
[see Eq. (123)]. We see that the drift is directed along the 
vorticity gradient. Assuming that the background vortic- 
ity u/(r) is positive and decreases monotonically with the 
radius, we find that the test vortex ascends the gradient 
if -To > and descends the gradient if r < 0. The drift 
velocity 77 = (V)drift can also be calculated from a linear 
response theory [22] . In the case of a unidirectional shear 
flow, combining Eqs. (21) and (23) of [55], it is found that 



(V) drift 



r. 



1 



f\S\ \ dw 
2.1^)1^^ 2 (83) 



where S(y) = (V)'(y) is the local shear. For t — ► +00, 
we obtain an expression equivalent to Eq. (|8ip for a cylin- 
drical shear flow. These expressions for the drift velocity 
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are equivalent to those obtained by Schecter & Dubin [32] 
using the Eulcr equation. These authors showed that ex- 
pressions (|8ip and (|83|) are only valid for retrograde vor- 
tices (Jo > if ijj is positive and decreasing). The linear 
response theory is not correct for prograde vortices (on 
the basis of numerical simulations [32]) so that nonlinear 
effects must be considered. 

For a thermal bath (statistical equilibrium state), in- 
troducing Eq. ([6"T]) in Eq. l[79"]) . using the property of the 
J-function to replace Q' by 12, and using Eq. ([6T]) again 
with r instead of r' , we find that the drift term (179]) takes 
the form 



r, = -D(5r 



dtp' 
dr 



(84) 



where D is given by Eq. (|78|) in the general case and by 
Eq. ([80]) when the profile of angular velocity is monotonic. 
In vectorial form, the drift can be written ij = —Dpr^Vij}'. 
It is perpendicular to the relative mean field velocity 
(V) = z x V-0'- In the analogy with the Brownian mo- 
tion [22|23j . the drift coefficient can be seen as a sort of 
"mobility" . It is related to the diffusion coefficient and to 
the temperature (which takes negative values in cases of 
interest) by an analogue of the Einstein relation [22] : 



i = D(3r . 



(85) 



We note that the drift velocity f2_3] of a test vortex is the 
counterpart of the Chandrasekhar dynamical friction [27j 
experienced by a star in a cluster. 



4.3 Connection with the usual form of the 
Fokker-Planck equation 

For an axisymmetric system, the Fokker-Planck equation 
is usually written in the form 



dP 
~dt 



2r dr 



d f((Arf 



dr V At 



-P 
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r dr 



rP 



(Ar) 
At 



where Ar denotes the increment of position of the 
test vortex in the radial direction. The second moment 
((Ar) 2 ) /2At represents the diffusion coefficient and the 
first moment (Ar)/At represents the drift velocity (the 
corresponding term in the Fokker-Planck equation can be 
viewed as an advection term). Comparing Eq. (|86p with 
Eq. ([711). we find that 



and 



D 



V 



1 ((Ar)' 

2 At 



(87) 



(Ar 



At 



dD 

dr 



The second expression shows that the drift velocity 
(Ar) I At is not exactly given by Eq. ([81]) but that there is 



an additional contribution dD/dr (sec [46]). Integrating 
Eq. ([75]) by parts, we have 



ri=— dr'tu(r')^-5(f2- f2')\n 

Ar J dr' 
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(89) 



Inserting the expressions ([75]) and ([59"]) in Eqs. ([57]) and 
([55]) . we find that the second (diffusion) and first (drift) 
moments of the position increment of the test vortex can 
be put in the form 



((Arf 
At 



7 /" +0 ° 

= '- / r'dr'5(f2 - Q') 
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(90) 
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V r dr 
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1 



(91) 



In Appendix [C] it is shown that these expressions can 
be obtained directly from the equations of motion of the 
point vortices in the large N limit. Therefore, starting 
from these expressions ([90]) and ([9T]) . inserting them in 
Eq. ([86]) . and repeating the above calculations in the other 
way round, we can derive the Fokker-Planck equation (IGOD 
for the relaxation of a test vortex in a bath. If we now ac- 
count for the fact that the distribution of the bath in Eq. 
(fBH]) is not stationary but slowly evolves under the effect 
of collisions in a self-consistent way, we obtain the integro- 
diffcrential equation ([24]) for the evolution of the system 
as a whole. This approach provides an alternative deriva- 
tion of the kinetic equations directly from the equations 
of motion (see Appendix [C| . 



(86) 4.4 Scaling of the collision term with N 



Let us consider the kinetic equation (|29p for a single 
species of point vortices. We introduce the total circu- 
lation r = Nj, the typical radius of the system R (fixed 
by the domain size or by the angular momentum) and the 
dynamical time t^ = R 2 / T. If we measure time in units 
of to and distances in units of R, the kinetic equation de- 
scribing the evolution of the iV-body system as a whole 
becomes 



duj 
~dt 



lid f +oa 
ANrdrJ y ' 



x In 



1 



, 1 dw 

r dr 



1 duj' 
r' dr' 



(92) 



This expression shows that the collision term scales like 
1/N. This implies that the kinetic equation ([9"2"]) cor- 
rectly describes the dynamics of the point vortex gas on 
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a timescalc ~ Ntjy (without h\N correction). We empha- 
size that this timescale does not generically represent the 
collisional relaxation time t re i ax towards statistical equi- 
librium because the kinetic equation (|92|1 does not relax 
towards the Boltzmann distribution in general. However, 
on this timescale, the system reaches a steady state with 
a monotonic profile of angular velocity (in generic situa- 
tions) due to distant collisions (resonances). Therefore, it 
is a sort of collisional relaxation time that we shall denote 



Nt, 



(93) 



where the asterix indicates that this is not, in general, the 
relaxation time t re i ax towards statistical eauilibriurj^l. 
The proper determination of the relaxation time towards 
the Boltzmann distribution (exactly), and the scaling of 



with N, remains unknown. 



Using the same normalization as before, the Fokker- 
Planck equation (|68|) describing the relaxation of a test 
vortex in a fixed distribution of N field vortices can be 
written 



dP _ 1 In N 1 d 
~dt ~ 8~r 



8P rflnH 



, (94) 



with a diffusion coefficient D(r) = uj(r)/\rS?'(r)\. This ex- 
pression shows that the typical timescale of the relaxation 
of the test vortex towards the bath distribution is 



t 



relax 



N 



to- 



rn 



Said differently, the relaxation time corresponds to the 
typical time needed by the vortex to diffuse on a length 
~ R. the domain size. Thus t re i ax ~ R 2 /D. Using Eq. 
(US}, we get t relax ~ R 2 \E\/(-y In Aw) ~ i? 2 /( 7 lniV) ~ 
NR 2 /(r\nN) ~ (N/ \nN)t D . We recall that the In N fac- 
tor arises from the logarithmic divergence of the diffusion 
coefficient when we make a bath approximation. There 
is no such term in Eqs. (f9"2"| and when we study 
the evolution of the system as a whole (this is different 
from the case of stellar systems where the In N term arises 
both in the gravitational Landau equation describing the 
evolution of the system as a whole and in the Kramers- 
Chandrasekhar equation describing the relaxation of a test 
particle [7]). In Sec. 15.21 we shall see that the spectrum of 
the eigenvalues of the Fokker-Planck equation has no 
gap when the diffusion coefficient decreases rapidly with 
the distance. This implies that the relaxation towards the 



12 There may be situations where the steady state reached 
by the system has an approximate Boltzmann distribution. 
Indeed, resonances have the tendency to "push" the system 
towards the Boltzmann distribution (since the Boltzmann en- 
tropy monotonically increases). Thus, if there are sufficient 
resonances during the evolution, the system can converge to a 
distribution close to the Boltzmann distribution on a timescale 
~ Ntn- This is however not always the case. For example, in 
the numerical simulation reported in Sec. 13.41 we have checked 
that the final distribution is relatively far from the Boltzmann 
distibution, even in the core. 



steady distribution (|71|) is not exponential but rather alge- 
braic. Therefore, Eq. (|95p is not an exponential relaxation 
time. However, it provides an estimate of the timescale on 
which "collisional" effects take place to drive the distri- 
bution of the test vortex towards the distribution of the 
bath. 



5 General study of the Fokker-Planck 
equation 

5.1 Relaxation of the tail of the distribution 



Equation (|68|) is a particular case of the general Fokker- 
Planck equation 



dP 
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(96) 



with d = 2, U(r) = — fjthi\u)\(r) where fx = /0/7, and 
D(r) is given by Eq. (|69|) . The stationary solution of this 
equation is 



P e (r) = Ae 



-U(r) 



(97) 



where A is a normalization constant. In an appropriate 
system of coordinates, the relaxation of the tail of the dis- 
tribution function has a front structure. This has been 
studied by Potapenko et al. [5SJ for Coulombian plas- 
mas and by Chavanis & Lemou fH] in the general case. 
Let us briefly recall the main results of the theory before 
considering its application to the case of point vortices. 
If we make the change of variables dx/dr = 1 / y/D(r) 
and introduce the function u(r, t) = P(r,t)/P e (r), we can 
rewrite the Fokker-Planck equation (|96jl in the form of an 
advection-diffusion equation 



at 

with a velocity field 



- + V[rix)] 0x-^ 



V(r) = -y/D(r) 



dr 



In 



r rf-l e -£7(r )£) l/2 (r) 



(98) 



(99) 



The evolution of the position of the front rj(t) is deter- 
mined by the equation 



^ = ^y)V(r f ). 



(100) 



Introducing z = x — Xf(t), u(x,t) = 4>(z,t) and r = 2t, 
the profile of the front for sufficiently large times is given 

by 



(z,t) = $ 



X(r) 



where 



+00 

<P(x) = -j= I e~ y "dy, 



(101) 



(102) 
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is connected to the error function by <P(x) = |[1 — erf(x)] 
and x( T ) is the function 



X 2 (t) = 2 J e W- H <r')W, (103) 
where H(t) is a primitive of h(r) = g(r/2) with 

g(t) = V'(r f (t))^D(r f (t)). (104) 

When g(t) = — l/(2t), a case that often occurs, we have 
x(t) = r 1 / 2 . The derivation of these formulae is given in 
|14j . In addition, it is shown that the parameter controlling 
the validity of the approximations made in the theory is 



e(t) 



4 V'(r f )/V"(r f ) 



(105) 



The theory becomes asymptotically exact if e(t) — > for 
t — > +oo. If this condition is not fulfilled, the theory can 
nevertheless provide a good description of the front struc- 
ture provided that e is sufficiently small (see explicit ex- 
amples in [14]). 



5.2 Correlation functions 

In this section, we determine the asymptotic behavior 
of the temporal correlation functions (A(Q)A(t)} of the 
Fokkcr-Planck equation ([96)) for t — > +oo. When the spec- 
trum of the linearized Fokker-Planck equation ([96]) has a 
gap (j, at the spectral value 0, the correlation functions 
tend to zero exponentially rapidly as e _Alt . However, when 
the spectrum has no gap in [0, +oo[, the decay can be 
slower (algebraic). This is the case in particular when 
the diffusion coefficient D(r) decreases rapidly with the 
distance (see below). This problem has been considered 
by different authors in various contexts [59-62,11] and we 
shall mainly rely upon their results. In particular, the sit- 
uation that we investigate here is closely related to the 
one considered by Bouchet & Dauxois [11] for the HMF 
model. 

In the following, we briefly recall the general theory 
developed in the Appendix B of Marksteiner et al. [59] 
(see also [60-62,11]) to compute the asymptotic behaviour 
of a correlation function associated to a Fokker-Planck 
equation with logarithmic potential and derive new re- 
sults that will be useful in the sequel. If we make the 
change of variables dx/dr = 1/ \J D{r) and f(x,t) = 
■J D(r)P(r, t)5 ( ir d_1 (where Sd denotes the surface of a 
d-dimensional sphere of unit radius), we can rewrite the 
Fokker-Planck equation ([9"B"]) in the form of a ID Fokker- 
Planck equation with a constant diffusion coefficient 



dt dx \ dx dx 



and a potential 
<2>(x) 



In 



r d-l e -t/(r)£>l/2 (r) 



(106) 



(107) 



We note that $'(x) = V^[r(x)] where V(r) is the velocity 
field introduced in Eq. I|99p. 

We recall that, in our context, the Fokker-Planck equa- 
tions are written in radial coordinates r > or x > 
and that the current must vanish at the origin. Therefore 
equations ([96]) or (|106[) are subject to the following radial 
boundary condition at r = or x — 0: 



|!E=0 = 0. 



(108) 



This is equivalent to say that the resulting solution can be 
extended to an even solution in ] — oo, +oo [. In particular 
the total mass on r, x > is conserved and a normalization 
condition can be assumed 



/(x, t)dx = 1. 



(109) 



Instead, the analysis developed by Marksteiner et al. [H] 
concerns the case of odd functions. This means that they 
consider a similar Fokkcr-Planck equation in which the 
boundary condition at r = is: 



P(t, r = 0) = 0, or f(t, x = 0) = 0. 



(110) 



In this case, the mass on r, x > is not conserved and the 
distribution function f(x,t) goes to for large time. 

We now give the main steps of the derivation of the 
asymptotic behavior of the temporal correlation functions 
in the natural context of the present study (this means 
that we incorporate the boundary condition Q108P with 
the normalization condition (|109[1 ) and compare the theo- 
retical results with numerical simulations. The stationary 
solution of Eqs. (|106|) - (|108|) can be written 



(111) 



where Z = J +o ° e~^ x ^dx is the normalization constant. 
We shall assume that <P(x) ~ a In a; for x — * +00. This 
is indeed a situation that occurs in the kinetic theory of 
point vortices (see below). The stationary solution then 
behaves like f e ~ x~ a for x — > +00 and it is normalizable 
provided that a > 1. Let W(x, t; xq, 0) be the solution of 
Eq. (|106|) with the initial condition f(x,0) = 5(x — xq)- 
For any function A(x), the temporal correlation function 
is defined by 



(A(O)A(t)) = J A{x o )A(x)W(x,t;x o ,0)f e (x )dx dx. 

(112) 

This can be rewritten 



+00 



(A(O)A(t)) = / A(x)f A (x,t)dx, (113) 



where /a(x, t) is the solution of Eq. (|106[) with initial con- 
dition Ja(x,0) = A(x)f e (x). Note that its total "mass" 
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Fig. 3. Evolution of the correlation function (|123[1 obtained 
by solving the Fokker-Planck equation (|106|) with a logarithmic 




potential ${x) = alnx for x > 1 and $(x) = ^a(x 2 — 1) for x < 



1. We have used the boundary condition (|108|l corresponding 
to even functions f(x,t). The correlation function is defined 
by Eqs. ([m]) and (fTT3)l with A(» = In as (8 = 1) for x > 1 
and -A(a;) = for x < 1. We have considered different values 
of a = 2, 3, 4, 5 (for clarity, the curves have been shifted to the 
same origin). To get rid of the constant (A), we have computed 
the time derivative C'(t) of the correlation function (this is a 
good strategy otherwise it is not clear whether we should use 
the exact or the numerical value of (A) which differ due to finite 
resolution) . The slope obtained numerically clearly depends on 
a and is found to be in good agreement with the theoretical 
prediction (|123[) . up to logarithmic corrections. 



is / + °° Ja{x, 0)dx = (A) where (A) is the average value 
of A(x) at equilibrium. This implies that /a{x, +oo) = 
(A)f e (x) since the mass is conserved. 

Setting -0 = /e*' 2 , we can transform Eq. (| 106[) into a 
Schrodingcr equation with imaginary time 



dip 
~dt 



d 2 il> 
dx 2 



(114) 



with a potential V s (x) = j^') 2 — \®" ■ Looking for solu- 
tions of the form ip ~ e~ Ekt ip k (x) with E k = k 2 (k > 0), 
we obtain the eigenvalue equation 

V s (x)ipk = Ekipk, 



d 2 jj k 
dx 2 



(115) 



with V s (x) ~ j/x 2 for x — ► +oo where 7 = a(a + 2)/4. 
Then, the general solution of Eq. (| 1 06[) with initial condi- 
tion /a(x,0) = A(x)f e (x) can be written 



f A (x,t)=e-*W* 



a ip (x) 



+00 



a(k)%p k (x)e 



-Ekt 



integrating on x and using the orthonormalization condi- 
tions, we obtain after simple calculations 

«(*0 = ~ / A(x)e-*W*it> k (x)dx, (117) 
z Jo 



dk 



(116) 



We assume that the functions ipk(x) are properly or- 
thonormalized, such that J ipk(x)ipk' (x)dx = 8(k — k') and 
/ ipo(x)ipk(x)dx = Skfi- Taking t = in Eq. pi6p . com- 
paring with the initial condition, multiplying by [x) , 



ao = - 



A^e-^^M^dx- 



(118) 



On the other hand, considering the limit t — > +00, we 
find that i/> (x) cx j e {x)e^l 2 oc e"*^'/ 2 . Therefore, the 
normalized ground state is 



ip (x) 



1 



-<S(x)/2 



(119) 



Inserting this expression in Eq. (|118|) . we find that 

a = J= ^ + °° A(x)/ e (x)^ = (120) 

Inserting Eq. pT6|) in Eq. pT3)l and using Eqs. (fTT7|) and 
(|120p . we finally obtain 



C(t) = (A(O)A(t)) (A)' 



Z I a{k) z e 
'0 



2„-E k t 



dk. 



(121) 



To obtain the large time behavior of the correlation 
function, we need to determine the equivalent of a(k) for 
k — > 0. To that purpose, we need the form of the eigenfunc- 
tions ipk(x) for k — > 0. In the case where V s (x) ~ 7/x 2 for 
x — > +00, corresponding to the logarithmic Fokker-Planck 
equation, it is possible to solve the eigenvalue equation 
(|115[) by match asymptotics |59lllj and obtain ipk{x) for 
k — > 0. Then, we can determine how the correlation func- 
tion ()121|) decreases for f — > +00 depending on the behav- 
ior of the function A(x) for x — > +00. Using the results of 
Marksteiner et al. [59], we find that when A(x) ~ x n the 
correlation function decreases algebraically like 



C(t) - t 



1 



(122) 



and 



This result has been given previously by Lutz [62| 
Bouchet & Dauxois [TT]. When A(x) ~ (lnx) 1 / 5 , using 
the results of Marksteiner et al. [59] , we find that the cor- 
relation function decreases like 



C{t) 



(lni) 



2/(5 



t~ 



(123) 



This expression has been checked numerically for different 
values of a (see Fig. [3]). For a = 3, we recover the result 
C(t) - (\ntf/ s /t given by Bouchet & Dauxois [EE] which 
turns out to be valid only for this particular value of a. 

The above results can be easily adapted to the case 
where the boundary condition is given by Eq. (|110|) in- 
stead of Eq. (|108[) . In that case, f e (x) is not a stationary 
solution of Eqs. ()106|) - (|110|) because it does not satisfy 
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Fig. 4. Same as Fig. [3] except that we solve the Fokker-Planck 
equation with the boundary condition (|110p corresponding to 
odd functions f(x,t). We have considered different values of 
a = 2,3, 4, 5, 6. In the present case (A) = so that we can 
directly plot C(t). The slope obtained numerically clearly de- 
pends on a and is found to be in good agreement with the 
theoretical prediction (|123[) . This is the same value as in Fig. 
[3] but we note that the scaling regime is reached much more 
rapidly. 



the boundary condition at r = x = 0. The stationary so- 
lution of Eqs. (fTO5] ) -([TTC3 ]) is f(x) = 0. This implies that 
ipo(x) = and ao = 0. As a result, Eq. (|121[) remains 
valid provided that we take (A) = 0. Furthermore, since 
the asymptotic behaviors of Eqs. (|122p and Q123[) obtained 
from Eq. (|121| only depend on the large x behavior of 
the functions ipk(x) and not on the boundary condition at 
x = 0, these results are valid both for the boundary condi- 
tions (|110p and ()108|) . This has been checked numerically 
in Figs. [3] and 21 Interestingly, these numerical simulations 
show that the scaling regime is reached more rapidly when 
the boundary condition is given by Eq. (|110j) instead of 
Eq. (fing]) . 



6 General asymptotic results 

We shall now apply the preceding results to the case of 
point vortices described by the Fokker-Planck equation 
([68]) . We first give general results obtained by considering 
the asymptotic behavior of the diffusion coefficient. 



6.1 Asymptotic behavior of the diffusion coefficient 

The angular velocity is determined from the vorticity by 
solving the differential equation 



uj = - — (r z Q). 
r dr 



(124) 



If Hi is a particular solution of this equation, the general 
solution is 0(r) = fi\(r) + C/r 2 where C is an arbitrary 
constant. One could a priori consider profiles of angular 
velocity that are divergent for r — > 0. However, this sit- 
uation is not physical and in general, the constant C is 



determined in order to avoid the singularity at r — 0. In 
that case, the profile of angular velocity is given by 



1 f r 

fi{r) = — / uj(s)s ds. 
r Jo 



Asymptotically, we have for r — > +oo: 



27rr 2 ' 



(125) 



(126) 



where r is the circulation. Then, using Eq. ()69|) . we obtain 
the asymptotic behavior of the diffusion coefficient 



D 



njlnA 2 2 

— — r u){r) oc r 0J{r). 



(127) 



More generally, combining Eqs. (|69[) and (|T24j) , we find 
that the diffusion coefficient can be expressed in terms of 
the angular velocity by 



£>(r) 



7, Ar 2 ny 



(128) 



Alternatively, for a given expression of the diffusion coef- 
ficient D(r), the above equation is a first order differential 
equation for (2(r). It can be easily integrated, leading to 



r2 ( r)=exp{-/ 



2dr 



ieD(r) I 
7 In A -I 



(129) 



where e = ±1 depending on the sign of the shear (e = +1 
if Q' > and e = -1 if O' < 0). We note that if D(r) 
for r — > +oo, which is the situation of physical interest, 
we get Q(r) ~ C/r 2 for r — > +oo. 



6.2 Power law 

We first assume that uj oc r~ a for r — > +oo. The circula- 
tion is finite if a > 2 and the angular momentum is finite 
if a > 4. The diffusion coefficient and the potential behave 
like 



D ~ D r 



,-(a-2) 



U(r) ~ fiahir. 



(130) 



In the following, we shall take Dq = 1 which can always 
be achieved by rescaling the time appropriately The sta- 
tionary solution (]9Tj) of the Fokker-Planck equation (|M|) 
decreases like P e ~ r _QM . It is normalizable provided 
that a/i > 2 and the variance (r 2 ) is finite provided that 
a.[i > 4. 

Let us first study the properties of the front struc- 
ture, using the results of Sec. 15.11 The change of vari- 
ables dx/dr = 1/ y/D(r) leads to r ~ (ax/2) 2 /" for 
r — > +oo. On the other hand, the velocity field (|9"9"]l 
behaves like V(r) - [(2fx + l)a/2 - 2]r" Q / 2 . Accord- 
ing to Eq. (|100p . the position of the front evolves like 
r f (t) ~ [a{{2fx + l)a/2 - 2)t] 1 / a . On the other hand, the 
function defined by Eq. (|104j) behaves like g(t) ~ — l/(2t) 
for t — > +oo. Therefore, according to Eq. (jlOip . we find 
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that the profile of the front characterizing the tail of the we find that the profile of the front characterizing the tail 



distribution for large times is given by 
u(r,t) ~ $ 



r a/2 



(2p+l)q 



2]t 



(2i)Va 



(131) 



For t — » +00, the parameter e(i) controlling the validity 
of the theory tends to 




4tt 



a((2/z + l)a - 4) 



(132) 



Since e(t) does not vanish for t — > +00, the theory is only 
marginally applicable. However, it can still provide a good 
description of the front structure when e is sufficiently 
small [H]. 

Let us now study the temporal auto-correlation func- 
tion of the position (r(O)r(i)}, using the results of Sec. 15.21 
The potential defined by Eq. p07p behaves like <P(x) ~ 
[((2fi + l)a — 4) /a] hxx. On the other hand, in terms of 
the variable x, we need to determine the correlation of 
A(x) = r(x) ~ x 2 / a . According to Eq. (TTJ5]), we find that 



(r(O)r(t)) 



_( ii o-4) 



(133) 



2 pre" 



of the distribution for large times is given by 
u(r,t) = <P 



X6 



(2t)Va 



(137) 



For t — > +00, the parameter e(t) controlling the validity 
of the theory tends to 



1 

e = — 



4~ 



8 V 2/i + 1 



(138) 



Since e(i) does not vanish for t — > +00, the theory is only 
marginally applicable. However, we shall see in Sec. 17.11 
that the theory still provides a good description of the 
front structure. 

Let us now study the temporal auto-correlation func- 
tion (r(O)r(i)) of the position, using the results of Sec. 
EH Since r ~ ((2/A) lux) 1 / 6 , the potential defined by 
Eq. <|107|) behaves like $(x) ~ (2/i + l)lnx. In terms of 
the variable x, we need to determine the correlation of 
A(x) = r(x) ~ (lux) 1 / 6 . Using Eq. (TT23"1) , we find that 



(r(O)r(t)} - (r) 



(Int) 



2/6 



(139) 



6.3 Stretched exponential 

We shall now assume that u ~ Ae~ Xr for r — > +00. The 
Gaussian case corresponds to <5 = 2 and the exponential 
case to 5 = 1. The diffusion coefficient and the potential 
behave like 



D 



D r 2 e~ 



U{r) - fiXr s 



(134) 



As before, we shall take Do = 1 without loss of generality. 
The stationary solution ([97]) of the Fokker-Planck equa- 
tion decreases like P e ~ e~ XfJ,r and all the moments 
of r are well-defined. 

Let us first study the properties of the front structure, 
using the results of Sec. 15.11 The relation between x and 
r, and the velocity field V(r) are asymptotically given by 



2 A a 

e 2r , 



XSr s 



rr( S (2/i + 1)A 5 _A r a 

k (r) ~ or e 2 . 



(135) 



According to Eq. (|100p . the position of the front evolves 
like 



e Xr s _ (2 M +1)A 2 



-25 



<S 2 t. 



(136) 



To leading order, we have rt(t) ~ (lni/A) 1 / 5 . On the other 
hand, the function defined by Eq. p04[) behaves like g(t) ~ 
— l/(2t) for t — > +oo. Therefore, according to Eq. (jlOip . 



6.4 Behaviour of position-correlations according to the 
considered vorticity profile 

We here qualitatively describe how the position- 
correlation function decreases depending on the consid- 
ered vorticity profile. According to Eq. (|106[) . the relax- 
ation towards equilibrium is controlled by the strengh of 
the drift term <P(x) ~ aolnx, i.e. by the exponent a$. 
Therefore, for large ao, the decorrelation should be fast, 
i.e. £ should be large, which is indeed the case according to 
Eqs. (|122p and (jT33J) . Now, if the vorticity profile decreases 
like cu(r) ~ r~ Q , we find that ao = 2/i+l — A/a. Therefore, 
ao is a monotonically increasing function of a tending to 
a = 2/i+l for a — > +00. On the other hand, the case 
w(r) ~ e~ Xr leads to a = 2/i+ 1 which is consistent with 
the case w(r) ~ r~ a with a — > +00. We conclude therefore 
that when the vorticity profile decreases rapidly (a large, 
or stretched exponential), the system decorrelates rapidly 
(ao large hence £ large) . The physical reason is not simple 
because the vorticity profile enters both in the diffusion 
coefficient and in the drift velocity in Eq. so that the 
decorrelation is a combined effect of these two terms. As 
discussed above, the situation is more easily interpreted 
from Eq. l|TP5Jt . 



7 Particular examples 

In this section, we treat explicit examples corresponding 
to typical distributions of the field vortices. 
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7.1 Gaussian vortex 

We first consider the case where the distribution of the 
field vortices corresponds to the statistical equilibrium 
state (thermal bath). As discussed in Sec.[2l it is obtained 
by maximizing the entropy of the point vortex gas 



S 



— In — dr, 

7 7 



(140) 



at fixed circulation, energy and angular momentum. This 
leads to the Boltzmann-Poisson equation 



(141) 



This maximization problem provides a condition of ther- 
modynamical stability. According to Appendix [XJ the 
functional (|140|) can also be interpreted as a particular 
H- function (or a Casimir functional). In this sense, a vor- 
ticity profile w(r) that maximizes S at fixed r, E and L 
is a nonlinearly dynamically stable stationary solution of 
the 2D Eulcr equation. We conclude that for the single 
species point vortex gas, thermo dynamical stability im- 
plies nonlinear dynamical stability. 

In this section, we consider the case where (3 — ► and 
Ql — > +oo in such a way that the product A = 7/3.!?l/2 
remains finite. The vorticity profile is then given by 

u = u e- Xr \ (142) 

This will be called the gaussian vortex (note that this is a 
particular case of Sec. 16. 31 with 6 = 2). The Fokker- Planck 
equation describing the relaxation of a test vortex in a 
gaussian bath is 



dP 
~dt 



ld_ 

r dr 



BP 

rD{r) ( — + 2\fiPr 



The stationary solution is 

P e (r) = 



V„-A M r 2 



(143) 



(144) 



It has been normalized such that J P e dr = 1. Using Eq. 
(|125p . the profile of angular velocity of the bath is 



2Ar 2 



(145) 



Therefore, according to Eq. (|69[). the diffusion coefficient 
is given by 



D{r) = 



D Xr 2 



i Xr2 - Ar 2 - 1 ' 



(146) 



with Dq = (1/4)7 In yl ( m the following, we shall take 
£>o A = 1 by rescaling the time appropriately). We note 
that the shear vanishes at r — leading to a divergence 
of the diffusion coefficient like D(r) ~ r~ 2 . This indicates 
a failure of the kinetic theory for r — > 0. For r — ► +00, 



The assumptions made in the kinetic theory assume a rela- 
tively strong shear so that the equations of motion for the point 
vortices are dominated by the mean field trajectories fTTJ). In 
the absence of shear, the expression of the diffusion coefficient 
is different as discussed in [2512] and [34]. Thus, the expression 
(|146|) is only valid at large distances r> 1. 



on the other hand, we have D(r) oc r 2 e~ Xr . Since we 
are particularly interested by the tail of the distribution, 
we shall extend this expression of the diffusion coefficient 
for all r, thereby circumventing the problems arising for 
r — > (this rcgularization shall not change the asymptotic 
results for r — > +00). 

Let us first discuss the front structure of the distribu- 
tion that is solution of the Fokkcr-Planck equation Q143P . 
According to Eq. (|136[) , the evolution of the front position 
is given by 



2(2/z+ l)X 2 t. 



(147) 



For very large times, it can be approximated by ~ 
(lni/A) 1 / 2 so that the evolution is very slow (logarithmic). 
On the other hand, according to Eq. (|137|) . the front profile 
is given by 



u(r, t) = <P 



1 ^ -V2(2/i+l)A^ 
A 



(2t)V2 



(148) 



These theoretical results are compared in Figs. [5] and [6] 
with direct numerical simulations of the Fokkcr-Planck 
equation (|143|) . For /i = 1, which is the value considered in 
the numerical simulations, the parameter (|138j) controlling 
the validity of the theory is e = i(47r/3) 1/2 ~ 0.256.... 
Since e(t) docs not tend to zero, the agreement with theory 
is not perfect, but the relatively small value of e explains 
why the theory gives however a fair description of the 
numerical results. In the present case, wc find that the 
front evolves less rapidly than predicted by the theory 
(assuming e — ► 0) but this tendency is not general. In 
P3] , we have found other examples where the front evolves 
more rapidly than predicted by the theory. Understanding 
this difference requires to develop a theory at order 0(e) 
in an expansion of the solutions of the front equation in 
powers of e — > 0. 

Finally, according to Eq. (|139[) . the temporal correla- 
tion function of the position of the test vortex is predicted 
to decrease like 



C(t) = (r(O)r(i)) - (rf 



lnt 

IF' 



(149) 



To check this prediction, the function (r(O)r(i)) is com- 
puted numerically from the expression 

(r(O)r(t)} = / r o rIT(r,t;ro,0)P e (r o )dr o dr 



J rP R (r,t)dr, (150) 



where Pu(r,t) = Pu(r,t) is the solution of the Fokker- 
Planck equation (|143p with initial condition Pr(t, 0) = 
rP e (r). On the other hand, (r(t)) is defined by 

(r(t)) = J rW(r,t;r 07 0)P e (ro)dr dr = J rP E {r,t)dr 7 

(151) 
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Fig. 5. Evolution of the front position r/(t) for the Fokker- 
Planck equation (|143p with a diffusion coefficient D(r) = 
r 2 e _Ar (we have taken fx = A = 1). The details of the pro- 
cedure are given in [14]. The result obtained numerically is 
in good agreement with the theoretical prediction (|147|1 . The 
agreement is, however, not perfect because the asymptotic 
value of e ~ 0.256 is not exactly zero. As a result, the front 
position (solid line) evolves less rapidly than predicted theo- 
retically (dashed line). 




r 



Fig. 6. Evolution of the front profile u(r, t) for the Fokker- 
Planck equation (|143[) with a diffusion coefficient D(r) — 
r 2 e~ Ar (we have taken fj, = A = 1). The solid lines corre- 
spond to the numerical simulation and the dashed lines to the 
theoretical prediction. The initial condition corresponds to the 
Heaviside function: P(r, 0) = l/(7va 2 ) if r < a and P(r, 0) = 
if r > a with a — 2. 



where Pe(r, t) = Pe(r, i) is the solution of the Fokkcr- 
Planck equation (|143j) with initial condition Ps(r, 0) = 
-P e (r). Since Pe(r, t) = P e (r) for all times (this is the 
stationary solution), we find that (r(t)) = (r) is the av- 
erage of r with the equilibrium distribution (|144p . Thus, 
(r) = ^(7r/A(u) 1//2 . We also note that / Padr = (r) so that 
Pr{ t i = ( r )Pe( r )- Therefore, according to Eq. p50p . 
(r(O)r(i)) — ► (r) 2 for t — ► +oo. This is why we need to 
subtract the term (r) 2 in the correlation function (|149[) . 
To avoid numerical errors in the evaluation of (r) (due to 
finite resolution), we have chosen to study the derivative 
of the correlation function which is predicted to behave 
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Fig. 7. Time evolution of the derivative of the correlation func- 
tion C(t). It is obtained from Eq. (|150|l by solving the Fokker- 
Planck equation (|143|l with the initial condition P(r, 0) = 

(l/ir)re~ r (we have taken A = fj, = 1). For large times, we 
find that the asymptotic slope of the curve is 2.0704..., in good 
agreement with the theoretical value 2 (up to logarithmic cor- 
rections) . 
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Fig. 8. This curve suggests that, for large times, the correction 
to the algebraic decay of the correlation function is logarithmic 
as expected from the theory. 

like 

C\t) - (152) 

The results of the numerical simulations are represented in 
Figs. [7] and [51 Figure [7] shows that the correlation function 
decays algebraically with the correct exponent predicted 
by the theory (up to logarithmic corrections considered in 
Fig.©. 

In conclusion, we find that the relaxation of the distri- 
bution function P(r, t) to the equilibrium state (|144[) is pe- 
culiar. The evolution of the front position is very slow (log- 
arithmic) and the temporal correlation function decreases 
algebraically (up to logarithmic corrections) . These results 
differ from the traditional exponential relaxation in the 
usual Brownian theory. They are due to the rapid (gaus- 
sian) decrease of the diffusion coefficient with r. The same 
behaviours have been found for the HMF model |ll|14j . 
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7.2 Isothermal vortex in a disk 



We now consider the case of an isothermal vortex ()141|) 
in a disk of radius R and assume that Ql = 0. This 
corresponds to a particular value of the angular momen- 
tum Lo(E,r) for a given energy E and circulation r. 
The equilibrium profile is then determined by solving the 
Boltzmann-Poisson equation 



ld_ 

r dr 



dtp 
dr 



Ae~ ^ = ujoe-^-^, (153) 



with ^'(0) = and i/j(R) = 0. This equation can be solved 
analytically, see e.g. [3912] . The vorticity profile is explic- 
itly given by 



Ar 



irR 2 (i] + A) 



ij+4 R 2 



(154) 



where r\ = f3jr/(2n) is a normalized temperature. From 
this distribution, we can determine the mean field en- 
ergy. The caloric curve (3(E) is given in [2]. There ex- 
ists solutions for any value of the normalized energy 
e = 2ttE/(N 2 r 2 ) > 0. The corresponding temperature 
satisfies rj > r] c = -4, i.e. (3 > (3 C = Stt/^T). For > 0, 
the vorticity increases with the distance and for f3 < 0, 
the vorticity decreases with the distance. For (3 = (3 C , we 
have u(r) = rS(r) so that e — > +oo. Typical vorticity 
profiles are shown in [2j depending on the value of the 
temperature. 

The Fokker-Planck equation describing the motion of a 
test vortex in a thermal bath specified by the distribution 
CH is 



dP 
~dt 



ld_ 

r dr 



rD(r) 



dP 



( d 



A[iPr\r 



(155) 



The profile of angular velocity corresponding to the 
isothermal vortex (|154p is 



Q = - 



IP 



^ + 4)1--^ 



(156) 



Therefore, according to Eq. (f6"9"| . the diffusion coefficient 
is given by 



D(r) 



AR 2 



\v\ 



(157) 



Again, this expression is valid only for r>l. Setting x = 
(l/2)(r/R) 2 , C = -7]/(T] + 4), Do =-y\nA(r, + 4)/fo| and 
r = Dot/AR 2 , we can rewrite the Fokker-Planck equation 

as 



dP 
57 



d_ 

dx 



dP 

dx 



r acp 

~ 1 + 2Cx 



(158) 



with < x < 1. This is of the form of Eq. pTJ7j|) with 
a logarithmic potential <P(x) = 2/iln(l + 2Cx). With the 
change of variable ip = Pe®! 2 = (1 + 2Cx)^P , it can be 



transformed into a Schrodinger equation with imaginary 
time 



dtp d 2 tp 
dr dx 2 



with a potential 



V s (x) = 



AC 



(1 + Cx) 2 ' 



(159) 



(160) 



where A = 4/i(l + fi). Since the domain is bounded, we 
cannot apply the results of Sees. 15. I land l5. 21 which assume 
that the tail of the distribution extends to infinity. There- 
fore, in the next section, we shall consider the isothermal 
vortex in an unbounded domain. 



7.3 Isothermal vortex in an infinite domain 

In an infinite domain, extending the calculations of [2], 
we find that the Boltzmann-Poisson equation (|153j) has 
solution only for the critical temperature 



13 =-- 

He — p! 



(161) 



and that the vorticity profile is given by the family of 
solutions 



u(r) 



w 



(1 + 2^1^)2 • 



(162) 



parameterized by the central vorticity ujq. Since u> ~ r~ 4 
at large distances, we note that the angular momentum di- 
verges logarithmically. From the Poisson equation ([7|) with 
the vorticity field (|162[) . we find that the streamfunction 
is given by 



ijj(r) 
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(163) 



where we have assumed that ip + 



for r 



in order to determine the Gauge constant. The energy of 
interaction E = \ J unpdr is finite and its value is 
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(164) 



This relation determines the central density uiq as a func- 
tion of E. The caloric curve is a straight line with con- 
stant temperature (3 = (3 C for all energies. We note that 
the kinetic energy E' = | J u 2 dr = | J(Wtp) 2 dr diverges 
logarithmically. It differs from the energy of interaction E 
because the boundary terms do not vanish in the present 
case. 

The Fokker-Planck equation describing the relaxation 
of a test vortex in a thermal bath specified by the distri- 
bution (fT62|) is 



dP 
~dt 



ld_ 

r dr 



rD(r) 




A^Pr 



(165) 
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The profile of angular velocity corresponding to the 
isothermal vortex ()162|) is given by 



Q{r) 



u 



2 1 



Therefore, the diffusion coefficient is 

i 

D{r) = ^-\nA^ 

47TCJQ ' 



(166) 



(167) 



Setting x = (iru>o/2r)r 2 , Do = 7r7myl and r = 
DoLOot I Ar , the Fokkcr-Planck equation can be rewritten 



dP d [dP 4/j.P 



dr dx \dx 1 + 2x J ' 



(168) 



with x G [0,oo [. This is of the form of Eq. (fTUS]) with 
a logarithmic potential •P(x) = 2/^ln(l + 2x). With the 
change of variable ip = Pe®' 2 = (1 + 2.x)''P, it can be 
transformed into a Schrodinger equation of the form ([1591) 
with a potential 



V s (x) 



A 



(169) 



where A = 4/x(l + /i). In the following, we set 
("fr/AnuJo) In A = 1 and r/iruo = 1 which can always be 
achieved by rescaling the time and the distances appropri- 
ately. Then, the results of Sec. \5\ can be directly applied 
with a — 4. Concerning the tail of the distribution, the 
position of the front is given by 



and the profile of the front is 



u{r,t) ~ 



(2i)Vs 



(170) 



(171) 



As discussed in Sec. 15.11 the theory is only marginally 
applicable. However, it can give a fair description of the 
front if the parameter e = (3/8)-\/7r/2/i is sufficiently small 
(for \i = 1, we have e ~ 0.47). On the other hand, the 
temporal correlation function decreases like 



(r(O)r(t)} - (r) 



t 



-(*»-!) 



(172) 



We note that we must impose /i > 1 otherwise the variance 
(r 2 ) of the equilibrium state diverges. 



7.4 Polytropic vortex 

As discussed in Sec. the distribution of the bath can be 
any stable stationary solution of the 2D Euler equation 
with a monotonic profile of angular velocity, not neces- 
sarily the statistical equilibrium state. A general criterion 
of nonlinear dynamical stability is given in Appendix [AJ 
As an explicit example, let us consider an i?-function (or 
Casimir integral) of the form 



1 



{uj q - u)dv. 



(173) 



A vorticity distribution which maximizes this functional 
at fixed circulation and energy (and angular momentum) 
is nonlinearly dynamically stable with respect to the 2D 
Euler equation. The variational problem Q19ip leads to a 
relationship between the vorticity and the stream function 
of the form 



/?(?-!) 



q 



(174) 



We shall call these vortices polytropic vortices |63j because 
they are analogous to stellar polytropes in astrophysics 
[20] . The structure and the stability of stellar polytropes 
in d dimensions has been studied in |64j . For d = 2, we can 
easily adapt these results to describe polytropic vortices. 
Note that the functional Q173P and the corresponding dis- 
tribution p74p arc similar to the so-called Tsallis entropies 
and Tsallis distributions introduced in non-standard ther- 
modynamics |65| . However, in the present context, they 
have a completely different interpretation [63,66-69]. They 
correspond to dynamical, not thcrmodynamical, equilib- 
rium states. In particular, the maximization of S q at fixed 
circulation and energy is a condition of nonlinear dynami- 
cal stability for the 2D Euler-Poisson system, not a condi- 
tion of generalized thermodynamical stability. The distri- 
butions (|174p form just a particular family of stationary 
solutions of the 2D Euler equation that can sometimes give 
a good fit of the quasi-stationary state (QSS) resulting 
from an incomplete violent relaxation [70156169] . This dy- 
namical interpretation of the functional (|173[) is different 
from the interpretation given by Boghosian [71] in terms 
of Tsallis generalized thermodynamics. For q — ► 1, we re- 
cover the isothermal vortex (|141|) as a special case. For 
q = 2, we have 



S 2 



(175) 



up to an additive constant. Therefore, the neg-enstrophy 
£2 = — -/~2 can be interpreted as a particular i7-function 
(or Casimir integral). In that case, the relationship (|174[) 
is linear 



1 



it) = /! 



(176) 



The structure and the stability of these linear vortices 
has been studied in detail in [72|70j . In these studies, 
the enstrophy was interpreted as an approximation of the 
Millcr-Robcrt-Sommcria mixing entropy [47-50] in a limit 
of strong mixing (or low energy). Therefore, the results of 
|72I70| were presented as results of thcrmodynamical sta- 
bility for the process of violent relaxation (in a particular 
limit of the theory) . Alternatively, if we view the enstrophy 
as an iJ-function (see Appendix |A"[) . the results of [72170] 
can also be interpreted as results of nonlinear dynamical 
stability with respect to the 2D Euler equation. 

For illustration, let us consider the case of a disk of 
radius R and take flj, = 0. We also assume that the vor- 
tex fills the whole domain and that ui(R) = ip(R) = 0. 
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This implies that = in Eq. (|176| . We then obtain the 
Hclmholtz equation 



A0 + k 2 tp = 0, 



(177) 



where k 2 = —/3/2 (we have assumed (3 < 0). The solution 
is V = AJ Q (ar/R) where a = kR = 2.40482... is the first 
zero of the Bessel function Jq(x). Setting w = a 2 A/R 2 , 
the vorticity profile is 



to = u)oJo(ar/R). 



(178) 



The Fokker-Planck equation describing the relaxation of 
a test vortex in a bath specified by the distribution (|178jl 



dP 
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r dr 



oTq Ji(ar / B) 
jR J (ar/R) 



(179) 



Using (xJiY = xJq, the profile of angular velocity corre- 
sponding to the Bessel distribution (|178| is given by 



Q = ^lj l{ar/R) , 



(180) 



violent relaxation of the 2D Eulcr-Poisson system in the 
first regime of the dynamics (see discussion in the Con- 
clusion) . According to Appendix [XJ a maximum of the 
entropy (|183|) at fixed energy and circulation (thermo- 
dynamical equilibrium state in the theory of violent re- 
laxation) is nonlinearly dynamically stable with respect 
to the 2D Euler-Poisson system. For K — > +oo, we re- 
cover the isothermal vortex (|14ip and for K — > 0, we ob- 
tain a vortex with uniform vorticity surrounded by irrota- 
tional flow (this is the analogous of the Fermi distribution 
in phase space). Considering now the limit f3 — > and 
S7l — * +oo in such a way that the product A = pOt /2 
remains finite, we obtain 
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(185) 



For K — > +oo, we recover the gaussian vortex (|142|) . More 
generally, the bath distribution behaves like oJ ~ e _Ar for 
r — > +oo so that the results of Sec. l7.1l applv to the present 
context. 



Therefore, according to Eq. (p9")l . the diffusion coefficient 
can be written 



D(r) 



7 ^Mar/R) 

4 n 2J 1 (ar/R) - <fJ (ar/R) ' 



(181) 



Considering now the limit (3 — *■ and i?£ — *■ +oo in such a 
way that the product A = /3Ql /2 remains finite, we obtain 



w(r) 



2 



(182) 



For g — > 1, we recover the gaussian vortex (|142[) . For q > 
1, the vorticity vanishes at a finite radius i?* (wc take 
w = for r > R*). For q < 1, the vorticity decreases like 
w ^ r -2/(i-< J r) Thjg algebraic behavior falls in the class of 
solutions studied in Sec. 



7.5 Fermi-Dirac vortex 

As a last example, we consider a typical prediction of the 
statistical mechanics of violent relaxation of the 2D Euler 
equation [47-50]. In the two-levels approximation, the sta- 
tistical equilibrium state is obtained by maximizing the 
mixing entropy 
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at fixed circulation, energy and angular momentum. This 
leads to a distribution of the form 
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(184) 



with K > 0. This will be called the Fermi-Dirac vortex. 
This distribution of the bath can emerge from a complete 



8 Conclusion 

In this paper, we have developed the kinetic theory of 
point vortices in two-dimensional hydrodynamics. The in- 
terest of this kinetic theory goes beyond the realm of fluid 
mechanics since the point vortex gas constitutes a sys- 
tem of particles with long-range interactions (like self- 
gravitating systems) which displays peculiar features in 
regard to statistical mechanics and kinetic theory [2]. 
Starting from the Liouville equation and using the projec- 
tion operator formalism, we have derived a kinetic equa- 
tion (|29|) that describes the evolution of the system as a 
whole. This kinetic equation conserves all the constraints 
of the point vortex dynamics (circulation r = 0, energy 
E = and angular momentum L = 0) and increases the 
Boltzmann entropy monotonically (S > 0). However, it 
does not necessarily relax towards the Boltzmann distri- 
bution of statistical equilibrium because this kinetic equa- 
tion admits an infinite number of stationary solutions and 
the evolution stops when there is no resonance anymore 
(this happens when the profile of angular velocity becomes 
monotonic). In that case, the vorticity profile can remain 
frozen in a QSS with a non-Boltzmannian distribution un- 
til more complex (e.g., three body) correlations come into 
play. This has been shown numerically in Sec. 13.41 We 
have next considered the relaxation of a test vortex in a 
steady bath of field vortices. This relaxation is described 
by a Fokker-Planck equation (|60|) with a space dependent 
diffusion coefficient. We have studied analytically and nu- 
merically some properties of this Fokker-Planck equation 
regarding the evolution of the front profile in the tail of the 
distribution function and the decay of the temporal corre- 
lation function of the position of the test vortex depending 
on the distribution of the bath. For a thermal bath (gaus- 
sian vortex), we have shown that the evolution of the front 
position is logarithmic, scaling like rj(t) oc (hit) 1 / 2 and 
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that the temporal correlation function of the position de- 
creases algebraically (up to logarithmic corrections), scal- 
ing like (Ar(0)Ar(t)) oc Int/t. This is at variance with 
the exponential decay of correlations for the usual Brow- 
nian motion (with constant diffusion coefficient and lin- 
ear drift) but this is similar with results obtained for the 
HMF model |ll|14j . We would like to conclude this paper 
by discussing the different regimes that take place in the 
dynamics of point vortices. This discussion is important 
to clarify certain points in order to answer some of the 
criticisms raised by Dubin |34] and discuss specifically the 
domains of application of the different kinetic theories. 

Basically, a system of point vortices is described by the 
iV-body Hamiltonian equations of motion ((T])-([2|). From 
these equations, we can derive the Klimontovich equation 
for the exact vorticity field Lo exact (r, t) = ^ 7;<5(r — r^(t)) 
expressed in terms of (^-functions (see Eq. (15) of [23] ) 
or the Liouville equation for the iV-body distribution 
P/v(ri, tat, t) (sec Eq. (59) of [23]). These exact equa- 
tions are equivalent to the A^-body dynamics and contain 
too much information to be of practical use. The object of 
the kinetic theory is to derive an equation for the smooth 
vorticity field u(r, t) = (J2i liH r ~ r i(t))} 01 the one-body 
distribution Pi(r,t) (see Sec. 13. ip . By rescaling the vari- 
ables appropriately, we have seen that the "collision term" 
could be evaluated in an expansion in powers of l/N for 
N — > +oo with fixed normalized energy e = E/(N 2 j 2 ) 
and fixed normalized temperature rj = 0jT. This un- 
usual thermodynamic limit corresponds to N — > +oo with 
7 ~ l/N and V ~ 1 (implying E ~ 1 and (3 ~ N). For 
the moment, only the term of order 0(1/ N) in the kinetic 
theory has been computed [see Eq. ([T4])]. This amounts to 
neglecting three-body (or higher) correlations among the 
vortices. These points are made clear in |46| . 

To leading order in N —> +oo, the evolution of the 
smooth vorticity field uj(r,t) representing the spatial dis- 
tribution of the point vortices is governed by the 2D Euler- 
Poisson system (see Eq. (125) of [23]) as for inviscid and 
incompressible 2D flows with continuous vorticity. This 
Vlasov (or mean field) limit amounts to neglecting all the 
correlations among the vortices. Physically, this regime 
is valid for t -C Ntjy and its domain of validity can be 
very long when N is large. The 2D Euler-Poisson system 
has a very complex dynamics. Starting from an unstable 
initial condition uj (r), the 2D Euler-Poisson system de- 
velops a complicated mixing process leading ultimately to 
a Quasi Stationary State (QSS), also called a mctaequi- 
librium state, which has the form of a vortex or a jet. 
This "collisionless relaxation" is similar to the process of 
violent relaxation in astrophysics based on the Vlasov- 
Poisson system |47|50|2| . It takes place on a timescale of 
the order of a few dynamical times to- The resulting QSS 
is a stable stationary solution of the 2D Eulcr equation on 
the coarse-grained scale uJ(r), i.e. if we locally average over 
the filaments. A statistical theory has been developed by 
Miller [48] and Robert & Sommeria [49] in order to predict 
the most probable equilibrium state resulting from a vio- 
lent relaxation. This theory is based on a hypothesis of er- 
godicity, assuming that the relaxation is complete, so that 



the system is expected to be in the most mixed state. This 
is similar to the statistical theory of Lyndcn-Bell [37] for 
collisionless stellar systems. A kinetic theory of the process 
of violent relaxation has been proposed by Robert & Som- 
meria [73] , and further developed in [50174175] , by using a 
phenomcnological Maximum Entropy Production Princi- 
ple (MEPP). They obtained a drift-diffusion equation for 
the coarse-grained vorticity uJ(r, t) with a time dependent 
temperature /3(t) that evolves in order to satisfy the con- 
servation of energy (see Eqs. (11)-(12) of [73]). In order 
to determine the domain of validity of this phenomeno- 
logical approach and determine the value of the diffusion 
coefficient (which is not given by the MEPP), Chavanis 
[76] attempted to develop a kinetic theory of violent re- 
laxation, starting directly from the 2D Euler equation and 
using a systematic procedure based on a coarse-graining 
and a linearization of the equation for the fluctuations. 
This quasilincar theory was developed in analogy with the 
quasilincar theory of the Vlasov-Poisson system in astro- 
physics [77178] . The approach of Chavanis [76] remains, 
however, unsatisfactory because it yields an equation [see 
Eq. (20)] which does not conserve energy. This could be 
cured as proposed in [23] [see Eq. (148)] by taking into ac- 
count resonances between the motion of the fluid particles 
but this procedure may not be completely justified since 
there is no real timescale separation between the decrease 
of correlations and the evolution of the coarse-grained vor- 
ticity field (unlike during the slow collisional relaxation of 
point vortices). As a result, memory (non-Mar kovian) ef- 
fects are expected to be important in the process of violent 
relaxation. Note that the more general equation (18) de- 
rived in [76] may conserve energy and provide a correct 
description of the coarse-grained dynamics [46] . 

If we now consider an axisymmctric distribution of 
point vortices that is stable with respect to the 2D Eu- 
ler equation and if we take into account terms of order 
1 /N in the kinetic theory, we describe the "collisional" re- 
laxation process. The kinetic equation governing the evo- 
lution of the system as a whole is given by Eq. ([29]) ob- 
tained by Chavanis [23] . or by the more general kinetic 
equation taking into account collective effects obtained 
by Dubin & O'Neil [51] . Physically, this regime is valid 
on a timescale Nto- For single species systems, the evo- 
lution of the smooth vorticity field is due to a condition 
of resonance f2(r, t) = f2(r', t) which can be satisfied only 
when the profile of angular velocity is non- monotonia As 
a result, the evolution stops when the profile of angular 
velocity becomes monotonic so that there is no resonance 
anymore. In that case, the system remains frozen in a long- 
lived QSS. This is a stationary solution of the 2D Eulcr 
equation which usually differs from the statistical equilib- 
rium state (Boltzmann distribution) predicted by Joyce & 
Montgomery [37] and Lundgren & Pointin [38]. The fact 
that the system can be frozen in a QSS state in the ab- 
sence of resonance was stressed in [33] . This implies either 
that the system will not reach statistical equilibrium at 
all or that the kinetic theory is incomplete. Indeed, the 
collision term calculated in the kinetic theory of Dubin & 
O'Neil [31] and Chavanis [23] corresponds to a term of or- 
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der 1 /N in a systematic expansion carried out in a proper 
thermodynamic limit N — > +00. In particular, it neglects 
non-trivial three-body correlations. Therefore, this theory 
describes the dynamics of the system on a timescale of 
order Nto where to is the dynamical time. If the col- 
lision term cancels out because of the absence of reso- 
nance, i.e. when the profile of angular velocity becomes 
monotonic, the relaxation towards statistical equilibrium 
(if any) will take place on a timescale larger than Nto 
so that next order terms in the 1/N expansion must be 
taken into account in the kinetic theory. These higher or- 
der terms have not yet been considered. Therefore, there is 
no certitude, up to date, that the point vortex gas relaxes 
towards statistical equilibrium for t — > +00. This could 
be checked in numerical simulations of the point vortex 
dynamics by solving the Hamilton equations (H}-©. A 
challenging problem is to determine the collisional relax- 
ation time t re i ax of the point vortex gas towards statisti- 
cal equilibrium. According to the present kinetic theory it 
should scale like t re i ax ~ N s to with 6 > 1 but the precise 
value of 8 is not yet known (see Sec. 14. 4[) . In numerical 
simulations, we must be careful not to confound the true 
statistical equilibrium state resulting from a slow colli- 
sional relaxation for large times (larger that Nto) with 
the QSS resulting from a violent collisionless relaxation 
taking place on a much shorter timescale ~ to- Indeed, 
there are many confusions in the literature between the 
Joyce-Montgomery [37] statistical equilibrium state (colli- 
sional relaxation) and the Miller-Robert-Sommeria [48 49 
statistical equilibrium state (collisionless relaxation). Al- 
though these statistical distributions are mathematically 
similar in certain limits (dilute limit), their physical ori- 
gin is completely different since they apply to different 
regimes of the point vortex dynamics. This implies in par- 
ticular that the two theories are not in conflict since they 
refer to completely different regimes. 

Finally, if we consider the relaxation of a test particle 
in a thermal bath of point vortices at statistical equilib- 
rium, we obtain a Fokkcr-Planck equation (|73|) which has 
the form of a drift-diffusion equation with a space de- 
pendent diffusion coefficient inversely proportional to the 
local shear created by the bath. This Fokker-Planck equa- 
tion derived by Chavanis [22123] is similar to the relax- 
ation equation proposed phcnomcnologically by Robert & 
Sommeria [73] by using the MEPP to describe the evolu- 
tion of the 2D Euler equation on the coarse-grained scale 
during the process of violent relaxation. However, this re- 
semblance is essentially coincidental since the two theo- 
ries describe completely different situations and different 
timcscalcs. Therefore, the kinetic theory of point vortices 
is rich and different regimes can be evidenced. For that 
reason, it is important to delimit the domains of appli- 
cation of each kinetic theory in order to avoid misunder- 
standings. Indeed, the mathematical structure of the ki- 
netic equations can be formally similar while describing 
completely different situations. When the conditions of 
validity of the different theories arc properly delineated, 



1 Some interesting work in that direction has been done re- 
cently in [57179] . 



there arc no contradiction between our results and those 
obtained by Dubin [34] . The works are complementary. 

In future works, it will be important to evaluate terms 
of higher order in the development in 1 /N of the kinetic 
equation so as to describe the evolution of the system on 
longer timescales. This implies the evaluation of two and 
three-body correlation functions, which is a formidable 
task. We must also distinguish the case where the cir- 
culations of the point vortices have the same sign or a 
different sign. When the circulations of the point vortices 
take positive and negative values, the point vortices can 
group themselves in dipoles or tripoles. As indicated in 
Sec. [5] this situation cannot be described by the mean 
field theory. The system will then have a very different 
dynamics than the one described here (in Sec. 13 . 2[) . For 
example, Sire & Chavanis [80] show that the merging be- 
tween vortices in 2D decaying turbulence is due essentially 
to collisions between dipoles and monopoles. The inclusion 
of dipoles and tripoles in the kinetic theory needs to take 
into account two and three-body correlation functions as 
discussed by Newton & Mezic [81]. In our approach, these 
pairs would form on a timescale larger than Nto so they 
have not been taken into account. We note, finally, that 
our approach implicitly assumes that there are no cor- 
relations initially so that the 1/N development is well- 
defined (see [46] for more details). If pairs (+,+), dipoles 
(+,—), tripoles (+,—,+) or more complex "structures" 
exist initially, then we must take into account high order 
correlations since the start and the kinetic theory will be 
different. 

Another extension of the kinetic theory is to consider 
flows that are not axisymmetric. In that case, the kinetic 
equation valid at the order 1/N is given by Eq. (fT4"]h It 
describes the evolution of the system as a whole on a 
timescale Nto- Using the timescale separation between 
the dynamical time to and the collisional relaxation time 
treiax (of order Nto or larger), it would be of interest to 
derive an "orbit-averaged" kinetic equation in terms of 
appropriate variables similar to the angle-action variables 
used in other contexts [20I18| . The resulting kinetic equa- 
tion would be more complicated than Eq. (p2"9"]) but the 
basic concepts remain the same. For non-axisymmetric 
flows, we can have a lot of resonances (see the expect- 
edly similar case in [TS]) and they can "push" the system 
towards a distribution close to the Boltzmann distribu- 
tion on a timescale ~ Nto (the exact Boltzmann distri- 
bution should be reached however on longer timescales 
due to higher order terms 1/N 2 , 1/N 3 ,... in the kinetic 
equation) . The same remark applies to the inhomogencous 
phase of the HMF model as discussed in [IS] . While there 
is no evolution at all on a timescale Nto in the homoge- 
neous phase of the HMF model (because there is no res- 
onance), many resonances appear in the inhomogeneous 
phase which can "push" the system towards a distribu- 
tion close to the Boltzmann distribution on a timescale 
Nto- These interesting problems should be given further 
consideration in future works. 
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A Nonlinear dynamical stability of stationary 
solutions of the 2D Euler equation 

We consider the 2D Euler-Poisson system 
du 



u • Vw = 0, u = z x V-0, uj = —Aip. 



(186) 



Any relation of the form uj = f(ip) determines a station- 
ary solution of the 2D Euler equation dui/dt = since 
Vw = f'(ip)Vip is perpendicular to u. This result is valid 
in an arbitrary domain. In a disk or in an infinite domain, 
we can also consider pseudo-stationary solutions of the 
form uj(r,8,t) = uj(r,8 — fliX) that describe a structure 
rotating uniformly with angular velocity fl^ (the system 
is stationary in the rotating frame). If we write the 2D 
Euler equation in the form 



duj 1 dip duj 1 dip doj 



= 0, 



(187) 



dt r d6 dr r dr d9 
and substitute for the preceding relation, we find that 

dip duj ( ' dip\ duo „ ,_„„s 

-m*- )de=°- (188) 



ij?x,r 2 is the relative stream- 



This equation is satisfied by any relation of the form 
uj = f(ip') where tp' = ip 
function. 

One important question concerns the dynamical sta- 
bility of a stationary solution of the 2D Euler equation. 
Let us introduce the functionals 



S= - J C(u)dr, 



(189) 



where C is any convex function, i.e. C" > 0. Such func- 
tionals form a particular class of Casimirs. They are also 
called generalized _ff- functions for the reasons explained in 
[82183] . It can be shown that the maximization problem 

Max{S[w] | E[uj] = E, r[uj] = r, L[ui\ = L}, (190) 

determines a steady solution of the 2D Euler equation with 
monotonic uj = f(ip') relationship that is nonlinearly dy- 
namically stablc^J. This condition of nonlinear dynamical 
stability has been proven by Ellis et al. [84) . It refines the 
Arnold's sufficient criteria of nonlinear dynamical stability 
and solves some apparent "paradoxes" observed in geo- 
physical fluid dynamics (some flows can be dynamically 
stable according to the criterion (|190j) although they do 
not satisfy the Arnold criteria) |84I85| . Note that the non- 
linear dynamical stability criterion (|190|1 has been proven 
in the absence of angular momentum constraint. When the 
angular momentum is conserved, the optimization prob- 
lem (|190j) determines a family of solutions that can be 



15 Note that this maximization problem is formally similar to 
a criterion of (generalized) thermodynamical stability in the 
microcanonical ensemble but it has a completely different phys- 
ical interpretation [63 68 . 



deduced from each other by a mere rotation. In that case, 
the notion of nonlinear dynamical stability must be recon- 
sidered to account for this property. 

Introducing Lagrange multipliers and writing the first 
order variations as 



we get 



8S - /35E - \pn h 5L - aSr = 0, (191) 



C'(uj) = -(3ip'-a, (192) 



where ip 1 is the relative stream function defined above. 
Since C is monotonic, this relation can be reversed to give 
uj = F(j3ip' + a) = f(ip') where F(x) = {C')- l {-x). 
Therefore / is monotonic. It is increasing at negative 
"temperatures" and decreasing at positive "tempera- 
tures" , since u'(ip') = -p3/C"{oj). When (3 -> and H L 
+oo in such a way that the product A = remains 
finite, the previous solution becomes uj = F(Xr 2 + a). 

In conclusion, the critical points of the optimization 
problem (|190|) determine steady states of the 2D Euler 
equation (in a rotating frame). Furthermore, maxima of S 
at fixed E, r (and L) are nonlinearly dynamically stable. 



B The calculation of the function M 

The velocity (by unit of circulation) created by point vor- 
tex 1 located in ri on point vortex located in r is 



r — Fi 

V(l -> 0) = — z x 



2tt 



r-ri 



1 

2^ 



z x Vlnlr — ri| 



(193) 



Therefore, the radial component V r = V • r of the velocity 
in the direction of the test vortex can be written 

^0) = -^|ln|r- ri |. (194) 

Introducing a polar system of coordinates (r, 9) to localize 
a point vortex and using 

|r - ri | = r\ + r 2 - 2rr x cos(<9 - 0r), (195) 

we obtain Eqs. (fH)j) and ([2U)l . The function M can be eval- 
uated from these expressions as done in |23I2| . However, 
it is more convenient to directly use the expansion 



ln|r-ri| =lnr> - ^ — - ( — ) 

^— ' 2 m V ? *> / 
m ^o 1 



i) 



so that 



1 x - vm I r< 
Airr ^ \m\ \ r> 

rn=£0 



|m| 



o ira{0-9 1 ) 



(196) 



(197) 
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Using the equations of motion (|17j) . we have 



V r (l->0,t-T) 



1 



im I r< 



4?rr Iml I r> ' 

m^O 1 



m 



im(9-8 1 -Anr) 

(198) 



Therefore, the function (|2"2")l takes the form 

1 



x E 



M = - 

ran / r< 
|mn| \r> 



+ 00 />27T 

J 16n 2 r 2 

\m\ + \n\ 

e im(e-ei) e in(e-6i- z!J7t) 



Integrating on the angles, we obtain 

1 / \ 2 l m l 

■ITT? 1 2 / ^ 



M 



dre 



imAQr 



Finally, the time integration yields 

2\m\ 



M =±y ( r <. 



m/0 v ^ 



8(mAf2), 



or equivalently 

- +oo / \ 2m 

^ 4,.2^( Z ^^) E ( r> J TO 



4t- 2 

This can be rewritten 
1 



M = ^5(AO) In 

Ar 2 



1-1^ 



(199) 



(200) 



(201) 



(202) 



(203) 



As pointed out in Sec. 13. 11 this expression has no sense 
for r and r\ such that fl{r) = fl{r\) and dfljdr(r\) = 0, 
because the Dirac function cannot be defined in this case 
even in the sense of distributions. However, a simple reg- 
ularization can be introduced to overcome this technical 
problem. Indeed, for the derivation of the expression M 
above, one could restrict the integration to a finite time t 
(this is in fact what the kinetic theory says at the start) 
and get an approximation M{t) of M as follows: 



M W = T7T~ J E 

m^O 

which is equivalent to 
1 



2|m| r +t 



dre mAnT , (204) 



M(t) 



E 



g 7rr 2 i — , y r ^ 

m>l 



^ 2m e imtAQ _ ^imtAQ 



iAfl 



(205) 



Using the complex logarithmic function, we get 
M(t) 



8nr 2 iAfl 









|-ln 


1 - 











In 



r < \ e -itAQ 



(206) 



which finally simplifies into 
1 1 



M(t) 



Airr 2 AQ 



arctan 



svaitAfl) 



1 - ( ^ ) cos(tZ\J7) 



(207) 



Note that this expression tends to expression ()203|) of M 
when t goes to infinity. 



C First and second moments of the radial 
increment 

In this Appendix, we calculate the first and second mo- 
ments (Ar) and ((Ar) 2 ) of the radial increment of the test 
vortex directly from the Hamiltonian equations of mo- 
tion |T])-([2]). We follow a procedure similar to that used 
by Valageas pjj in a different context. Since the calcula- 
tions are similar, we shall only give the main steps of the 
derivation. For simplicity, we assume that all the vortices 
have the same circulation 7. In order to separate the mean 
field dynamics from the discrete effects which give rise to 
the diffusion and the drift of point vortices, we write the 
Hamiltonian j2]) as 

H = 1 (H Q + H I ), (208) 
where we defined the mean field Hamiltonian Hq by 

N 



H =^^ (ri), 



and the interaction Hamiltonian Hj by 



-^E ln i^- r ii-E^) 

i^j i 



(209) 



(210) 



In Eq. (|209[) the mean field stream function is given by 

Mr) = ~ I w(r') In |r - v'\dv', (211) 

where u(r') is the mean field equilibrium vorticity. The 
factor e ut has been added for the computation of perturba- 
tive eigenmodes and we shall ultimately let uj — > + (this u 
should not be confused with the vorticity w(r)). Thus H 
describes the mean field dynamics whereas Hi describes 
the discrete effects which vanish in the limit N — > +00. 
Therefore, we consider Hj as a perturbation of H and we 
apply a perturbative analysis in powers of 1/N . 

Using the identity (|196j) . the interaction Hamiltonian 
Hi is given by 



4tt 



2\m\ \r> 



im(9-0') 



47T 



E E 2 k 

In r> + ^ ] J Lj(r')\nr > r l dr 1 



(212) 
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On the other hand, introducing a polar system of coordi- 
nates, the equations of motion read 

dn 18H d6i 8H 

7 dF = ^' in ^ = -^- (213) 

We write the trajectories {r(t),6(t)} as the perturbative 
expansions r = r^ + r ( - 1 - ) + r^ + ... where r' fe ' is formally 
of order k over Hi. At zeroth-order, we simply have 



dr 



(0) 



0, 



dO 



(0) 



^ (0) ), 



(214) 



dt dt 
which yields the mean field equilibrium orbits 

rf ] = constant, of ] = Q{rf ] )t + 0f ] (0). (215) 
At first order, we obtain 



dr 



(i) 



1 dH! 



dt r (o) dOi 



dO 



(i) 



// (oh (i) 1 d-Hj 



dt 



= Q'(rDr, 



Jo) dn ' 

i 

(216) 



where we can substitute the zeroth-order orbits in the 
r.h.s. A simple calculation yields 



(i) 



j_dx 

,(°) dOi ' 



,(°) dr* ; 



(217) 



with 



1 / r< \ |m| e im(e-e') 



47r ^ ^ 2|m| V. r-> / w + im(J? - J?') ' 



(218) 



At second order, we have 

dr„ (2) 1 ^ <9 2 i?/ (i) 



eft 



_(o) X! 



(o) 96 i i5r J 
'i 3 J 



E 



1 ^ d 2 ff 7 „m r, (1) 3ff 



(0) a6 /. a ^ j £ r (o)ja 961, 



(219) 



Considering a test vortex, and taking the average with 
the distribution w(r) over the orbits of other field vortices 
yields at order 1/N 



■ ■(2)\ 7 2u>t ® f ill i l\ 

r ') = —e — r dr uj(r ) 
' 8?r dr J v ' 



x V — I — 

2— 1 r 2 lr> 

m#0 



2|m| 



w 2 + m 2 (r?- fl>) 2 



8 



8tt / v 'dr 



E 

m^0 



1 /r 



r \r> 



2 1 m | 



; 2 +m 2 (l2-L7') 2 ' 



(220) 



Then, using lim^, ^quj / (lu 2 + x 2 ) = ivS(x), the limit u> — -> 
0+ gives Eq. 1[9"T]). On the other hand, from Eqs. l[2l~7|) - 
(I218p . the mean square radial displacement {{Ar) 2 ) reads 
at order 1/N: 

{{Ar)2) = 8^ e2Wtl / r ' dr>LU ^ 



><E 



2|m| 



i 



r>y lu 2 + m 2 (n - Q 1 ) 2 



x 1 + e 



luiAt n„uiAt 



2e u ' At cos[m(L? - fi')At]) . (221) 



The limit u> — > + now gives 



(Z\r) ! 



1 



m^0 



At / Anr 2 

2\m\ 



r'dr'uj(r') 



1 - cos[m(L? - Q')At] 
m 2 (n- f2') 2 At 



(222) 



Taking At — > +oo and using lim t _> +00 (l — coste) jtx 2 
ttS(x), we obtain Eq. (|50|) . 
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